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ABSTRACT 

A description of the Generation-I Engine Simulator, and the analog- 
computer circuit for reactor simulation, is presented. The results of the 
analyses of the computer -control system for achieving simulated nuclear 
engine characteristics are also incorporated. The test runs and test results 
are extensively described including an analysis of the data. 
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I. INTRODUCTION 

A. BACKGROUND 

During the NERVA Engine proposal effort of Aerojet m early 1961, 
Aerojet-General Corporation undertook a test program to gain operating exper- 
ience and experimental data pertaining to the operation of a liquid-hydrogen 
propellant-feed-system in a nuclear rocket engine. 

The principal parts of the experimental feed system used in this 
program consisted of a centrifugal-type, liquid-hydrogen pump driven by a 
two-stage, impulse turbine, both components possessed the performance char- 
acteristics required for a NERVA Engine employing the cold-bleed cycle. To 
demonstrate the feasibility of heat-capacity startup, turbine fluid was bled from 
the pump discharge line and passed through an aluminum mockup of a reflector 
segment. A highly pressurized source of gaseous hydrogen was available to 
supply power to the turbine when the thermal capacity of the aluminum heat 
exchanger was exhausted. 

In order to regulate both the power input to the turbine and the pump 
load in accordance with the startup conditions in a nuclear rocket engine, a small 
portable analog computer was set up to simulate factors of reactor kinetics and 
heat-transfer and fluid dynamics in both reflector and reactor. It was anticipated 
that pump discharge pressure would be fed into the computer, and computer out- 
puts transferred to the physical system by means of four servo-actuated, flow- 
control valves. 

The resulting system, called the Cold-Bleed Engine Simulation Test 
System (C~BEST)» was tested in March 1961 with several pump tests and two 
heat-capacity star't tests. Failure of two power supplies in the portable analog 
computer ruled odt the use of that equipment m that test series because of the 
lack of spares, and because the test schedule allowed insufficient time for re- 
pairs. Some resd-lts of C-BEST tests are included in this report since they are 
pertinent to the analysis of the Generation I Engine Simulator. 

After award to Aerojet of Contract SNP- 1 for design of the NERVA 
Engine, the C-BEST system was revised to form the Gener ation-I Engine Simu- 
lator System vdnc'h. was tested in September and October of 1961, to fulfill the 
contractual obligation of Task Item 10, Contract SNP-1. 

Major changes from the C-BEST configuration included insertion of 
a second-generation turbopump, placement of the heat exchanger m a vertical 
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position, addition of thermocouples in intermediate locations in the heat exchanger, 
provision for longer flow- straightening sections ahead of flow-measuring orifices, 
and improvements of the overspeed safety system. Based upon additional data 
pertaining to the dynamic performance of the reactor control loop, which had be- 
come available from analyses conducted for the NERVA Proposal effort, the 
original analog-computer circuit was revised. An additional power supply was 
added to eliminate the previously incurred overload condition, the two power 
supplies that had previously failed were repaired, and the portable computer was 
entirely reworked to reduce the possibility of ground loops. Furthermore, experi- 
mental frequency responses were run on the flow-control valve servo-actuators 
after the conclusion of the C-BEST tests. 

In consideration of the design revisions and the new information avail- 
able from both C-BEST and actuator tests, it was necessary to treat the 
Generation-I Engine Simulator as a completely new system so far as the analysis 
of system performance and control system design was concerned. 
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B. TEST PURPOSES AND OBJECTIVES 

One major purpose of the Generation I Engine Simulator was to 
demonstrate the feasibility of initiating turbopump "bootstrapping" by utilizing 
the thermal energy available in the propellant-feed system and, especially, in 
the reflector of the engine prior to the time radiation heating becomes significant 
(in the NERVA engine, all the propellant flows through the reflector). Lack of 
time and funding factors did not permit the construction of a full-scale simulated 
reflector. However, a properly designed reflector segment could be used to 
heat turbine-drive fluid because the thermal capacity of the segment would have 
the same proportionality to the total heat capacity of the reflector as turbine flow 
would have to the total flow through the reflector. 

The second major purpose of the test system was to gain data and 
operating experience with a full-scale, liquid-hydrogen propellant-feed system 
during a start transient; information that was directly related both to the design 
of the propellant feed system, and to the analysis and design of the NERVA 
engine control system. It was necessary to obtain knowledge concerning the 
effect of various pump inlet pressures on the start transient, and on the boot- 
strap capability of the pump; in addition, experimental data concerning the 
dynamic response of the propellant feed system to variations of pump load and 
turbine energy was required for the NERVA engine dynamic analysis. During 
engine startup, and after the engine reactor reaches a power level sufficient to 
vaporize the propellant in the reflector, the propellant feed system is subjected 
to an increasing pump load simultaneously with an increasing energy input to the 
turbine. In order to produce these two related events, and to investigate the 
effects on turbopump bootstrapping of a continuous rise in reactor power, a small 
portable analog computer was provided to simulate, in real-time, the effects of 
reactor heat transfer and fluid dynamics. Although the inadequacy of such a 
small computer was recognized, still, it was considered useful as an adjunct in 
achievement of the test objectives. Furthermore, since the NERVA program 
called for the design of a much more sophisticated engine simulator involving an 
analog computer, and since Aerojet had never before attempted partial system 
testing of an engine, use of this smaller computer was felt to be justified m 
acquiring practical knowledge and experience attending the problems of inte- 
grating physical components with an analog computer to achieve simulated system 
performance. 

Another purpose of the test system was to investigate the transient 
heat-transfer characteristics of the simulated reflector segment, as a supplement 
to a theoretical analysis of heat transfer to boiling liquid hydrogen at pressures 
below the critical point. The analysis was aimed at determining the mechanism 
of two-phase flow instability, if it existed, and to derive a set of equations de- 
scribing the boiling process that would be applicable to the analog computer study 
of the startup control problem. 
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C. DESCRIPTION OF TEST SYSTEM APPARATUS 

The general configuration and relative locations of components of the 
Generation-I Engine Simulator is shown in Figure I-C, 1. Propellant flow was 
initiated by pressurizing the LH 2 tank with gaseous hydrogen. A regulating 
system on the tank was provided to hold the tank pressure, as preselected for a 
given test, to within a specified tolerance of i0.5 psi. 

Fluid left the pump and passed through a vaned, flow-straightening 
section, and through a square-edged flow-metering orifice, 3.807 inches in 
diameter. Part of the flow entered the bleed line, and the rest went through a 
parallel section, one branch of which contained a 2. 35-in. -dia. , square-edged 
orifice ( A^fYo)’ and the other a 5-m. butterfly valve (the load valve) which could 
be positioned by a servo-controlled hydraulic actuator. The bypass orifice was 
provided to avoid any possibility of unloading the pump m the event of load-valve 
closure. Following the bypass section, the hydrogen entered a duct, 140-ft long, 
and was burned off at the end. 

The bleed-line flow passed first through a square-edged orifice, 

1 . 25 -in. -dia, (APbpO^’ an d then through the bleed-line valve. Downstream 
from the bleed-line valve, the fluid entered at the bottom of the heat-capacity 
heat exchanger placed vertically m the test stand. 

A check valve was initially installed at the heat-exchanger exit to 
prevent gaseous hydrogen (which entered the mixing section from the gas farm) 
from blowing back through the bleed line to the discharge line where it could, 
by vaporizing the liquid propellant, unload the pump. 

A mixing section was provided to allow thermal mixing of liquid and 
gaseous hydrogen, which enters, respectively, from the bleed line and from the 
line to the gas-farm. Since in the gas-farm, hydrogen gas was stored at 2600 
psig, gas pressure was reduced to, and maintained at, 800 psig, by three pres- 
sure regulators upstream of the APqH 2S fl° w -metering orifice (2. 15-m. 
diameter, square-edged). The gas-farm valve controlled the flow of gaseous 
hydrogen. 


The turbine supply line contained the turbine power control valve, a 
square-edged orifice, 2.25-m. -dia, and a turbine safety shutoff valve, m that 
order, and terminated at the turbine nozzle. 

The portable analog computer, which was used for reactor simulation, 
consisted of two consoles containing 18 dc chopper-stabilized operational ampli- 
fiers. Of the 18 amplifiers, four were used as error detectors, as indicated in 
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Figure I-C, 1, m the servo loops, and the remainder were used for reactor simu- 
lation. The Universal Sequence unit, part of the facility at Test Stand C-6, was 
used to transfer control signals between the computer and the test stand. The 
patchboard in the sequence unit allowed a flexible setup whereby individual servo 
loops could be easily checked out prior to test, and facilitated test operations by 
permitting control functions to be changed at the control console. 


Figure I-C,2 shows an instrumentation diagram of the system, while 
instrument ranges are given in Table I-C 7 1. Based upon studies made by the 
Instrumentation Dept, of the Test Division, Liquid Rocket Plant, AGC, it was 
concluded that the maximum expected error m pressure transducers is 1 0. 48% 
of the point, and includes errors due to transducer inaccuracies, as well as 
those in the data acquisition system. These same studies show that a maximum 
error of i 0. 44 F m temperatures was obtained with platinum-element resist- 
ance probes, however, the data from Generation-I tests indicates discrepancies 
of 1 to 4 m temperatures below -400 F. The platinum resistance probes were 
used to measure LH 2 temperatures, whereas copper-constantm thermocouples 
were used for higher temperature ranges, as indicated in Table I-C,l. The 
accuracy expected of the c-c thermocouples was i 1 F for temperatures down 
to approx. -300°F. 


The above mentioned studies indicate that turbine-type flow meters 
for service have maximum errors of i 0. 85%. However, the flow meter used m 
Generation-I engine simulator tests was volumetncally calibrated with water, 
therefore, as indicated above, the measured temperatures of LH 2 may have been 
in error by a few degrees. Thus, mass flow rates calculated according to flow 
meter readings would be 10 to 15% m error because of the error involved in 
picking LH 2 density from a T-S diagram. 
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TABLE I-C, 1 

INSTRUMENTATION FOR GENERATION- 1 ENGINE SIMULATOR* 


F unction 


fT 


fs 


fD 


Ap 


fDo 


fYo 

Ap 


fYo 


LLV 1 


LLVo 


BPo 

Ap 


BPo 


BFCVi 


Hei 


Heo 


gh 2 s 


Ap 


gh 2 s 


GH 2 CVi 


GH 2 CVo 


Range 

psig ^F lb / sec rpm 

0 - 100 

0 - 100 

0 - 1000 

0 - 200 

0 - 1000 

0 - 200 

0 - 1000 

0 - 1000 

0 - 1000 

0 - 500 

0 - 1000 

0 - 1000 

0 - 1000 

0 - 1000 

0 - 500 

0 - 1000 

0 - 1000 
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TABLE I-C, 1 (cont. ) 




Range 



F unction 

psig °F lb/ sec 

rpm Pot 


P 

ms 1 

0 - 1000 



P 

mso 

0 - 1000 



P 

TS 

0 - 1000 



A p TS 

0 - 500 



P Ti 

0 - 1000 



P Te 

0 - 100 



T fT 

-430 to -300 



T 

fFM 

-430 to -300 



T 

fS 

-430 to -300 



T 

£D 

-430 to -300 



T 

fYo 

-430 to -300 



T 

BPO 

-430 to -300 



T 

Hei 

-430 to -300 



T 

Heo 

-300 to + 100 



T 

gh 2 s 

0 to + 100 

H 



T GH 2 V0(0) 

0 to +100 



T 

msi 

-300 to +100 



T 

mso 

-300 to zero 
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Function 


Ts 


Ti 


Te 


MF 


LLV 


L 


TPCV 


BFCV 


GH 2 CV 

J FSV 


g 


w 


1-6 


LLYo 


fDe 


fDe 


N 


psig 


0 - 100 


Range 

°F lb/ sec 

-300 to zero 
-300 to zero 
-300 to zero 

0 - 100 


Amb to -400 
Amb to -400 
-300 to -430 


rpm Pot 


Pot 

Pot 

Pot 

Pot 

Pot 


-300 to -430 


0 - 30,000 
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TABLE I-C, 1 (cont. ) 


For analog computer operation 


F unction 

psig 

Range 

°F 

Remarks 

P fD 

0 - 1000 


Set gam to 0. 010 volts/psi 

AP fDo 

0 - 100 


Set gam to 0. 100 volts/psi 

P 

mso 

0 - 1000 


Set gain to 0.010 volts/psi 

T 

mso 

• 

+ 100 to -320 

Amplify with Kmtel Model 110 or 
111. Gain setting of 70. 


For instrument locations, see Figure I-C, 2. 
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LEGEND 
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@ TURBINE FLOW METER 
(5) FUEL SUCTION VALVE 
( 4 ) PUMP 
(?) TURBINE 
GH 


(6) LOAD VALVE 

( 7 ) BLEED-LINE VALVE 

(?) HEAT-CAPACITY HEAT EXCHANGER 
(?) CHECK VALVE 

@ MIXING SECTION 


(H) TURBINE CONTROL VALVE 
@ GAS-FARM CONTROL VALVE 
Q SERVO AMPLIFIER 
(14) TURBINE SHUT-OFF VALVE 


in so 


POSITION 

FEEDBACly 



I I 

PORTABLE ANALOG COMPUTER 


SIMULATOR RUN 
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Figure I-C,l 
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D. TEST SYSTEM OPERATION 

1 . Heat-Capacity Start Tests 


In the first series of heat-capacity start tests, the analog com- 
puter was inoperative. The gas-farm valve remained closed at all times during 
these tests, while the bleed-line valve remained open. The turbine -power control 
valve (TPCV) remained either completely open or at a predetermined position. 
Similarly, the load valve remained fixed at a predetermined position, such that 
the resistance to flow of the parallel pipe section corresponded, approximately, 
to the resistance of the nuclear engine. The flow control valve positions were 
maintained by the hydraulic servo-actuators which received position commands 
from the sequencing unit. As Figure I-C-, I illustrates, the set of relays oper- 
ated by the "simulator run" switch were not energized for heat- capacity start 
tests; consequently, they were in the normally closed position so that the fixed 
voltages V j to Vq were applied to the servo -actuator s , which were made posi- 
tional since the position feedback loops were also closed 

Prior to the start of a test, which was initiated by opening the 
fuel suction valve, the safety valve was opened to allow the suction line to be 
cooled by liquid hydrogen. Tank pressure could be maintained at any desired 
level to provide a net positive suction pressure (NPSP). A portion of the hydrogen 
flowing through the pump entered the bleed-line where it acquired additional heat 
m the heat-capacity heat exchanger on its way to the turbine. Sufficient thermal 
energy was imparted to the fluid by the pump, lines, and reflector segment (heat 
exchanger) to drive the turbine; However, as the suction line and pump were 
cooled to the temperature of_LH 2 , the load on the pump increased, while the 
energy available at the turbine decreased; consequently, turbopump speed rose 
to a peak, then decayed. 

The level of propellant in the tank was monitored by means of 
a capacitance probe and when the level in the tank dropped to 15%, shutdown 
(FS-2) was initiated. On FS-2, the second set of relays m the sequence unit were 
energized, a new set of input voltages, Vg - Vg, were applied to the servo ampli- 
fiers to position the valves for the shut down operation. 

2. Computer -Controlled Test 


The start of a computer controlled test was identical to the 
start of a heat-capacity start test. From the start at FS-1, the computer proc- 
essed the signals from the pressure and temperature transducers, but the flow 
control valves were held to pre-set positions by the sequencing unit, that is, the 
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simulator-run switch was not closed. Not until stable LH 2 appeared at the pump 
exit, assuring a load on the pump, was control of the test system transferred to 
the analog computer. Control was transferred at the control room console by 
closing the simulator-run switch which operated the relays m the sequencing unit. 
When these relays were energized, the signals from the error detecting amplifiers 
in the computer were applied directly to the servo-valve amplifiers. Simultan- 
eously, the positional servos were converted to velocity controllers by removal 
of the valve position f eedback signals. The reason for delaying computer take- 
over was to avoid the possibility of applying auxiliary power from the gas farm to 
the turbine while the pump was unloaded. 

The primary input to the computer was pump discharge pressure 
(Pjq), as a voltage. The computer circuit was set up to solve, in real time, the 
equations describing fluid dynamics, heat transfer, and neutron kinetics m the 
nozzle, reflector, and reactor. As the result of applying pump discharge pressure 
to the computer, computer voltages were generated which could be considered to 
represent real engine variables such as chamber pressure and temperature, total 
flow rate, and pressure and temperature at the inlet to the TPCV. 

3. Reactor -Simulator Control 


The Reactor Simulator included the entire control loop, one of 
the two primary engine control systems. The Engine Programmer (which, in the 
real engine, programs the reference inputs of chamber pressure and temperature 
for the engine control systems) was also simulated. On transferring control of the 
test system to the Reactor Simulator, the simulated Engine Programmer started 
a programmed rise m P c and T c references, and the variables in the Reactor 
Simulator followed this program. The other primary control loop, which took 
thrust-chamber pressure (P c ) as a reference, regulated P c through the TPCV, 
which controlled pump discharge pressure. Part of the latter loop was real, 
since it included the prototype TPCV and the turbopump, while part of it was 
simulated, since P c and P c _ re £ were generated in the computer, though the 
computed value of P c was a direct function of pump discharge pressure. 

Pressures and temperatures generated as voltages in the com- 
puter were references against which the actual values, appropriately transduced 
and scaled, were compared. Any errors that existed were reduced to zero by 
the four flow-control valve servo-actuators. One output of the computer was the 
analog solution for total weight flow rate (w r ), which reflected the resistance to 
flow caused by a rise in reactor power. The flow rate was compared to actual 
pump delivery (wj), as measured by the differential pressure drop across the 
orifice at the pump exit (Figure I-C, 1), and any error was reduced to zero by 
the servo actuated load valve, with the result that actual pump delivery closely 
followed computed flow rate. 
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The bleed-lme and gas-farm valves acted together to control 
pressure and temperature at the inlet to the TPCV. Control of temperature was 
accomplished by comparing actual fluid temperature (T mso ) with the analog 
solution for reflector outlet temperature resulting from heat transferred to the 
propellant in the reflector. Admission of warm gas to the test system, by the 
gas-farm valve, compensated for the lack of radiation heating in the test apparatus. 
Pressure control was obtained in a similar way by using the error between actual 
pressure (P mso ) and computed pressure (P ro ) to position the bleed-line valve. 

These three valves (TPCV, and gas-farm and bleed-lme valves) 
were necessary to transfer computer outputs to the test system, that is, to load 
the pump and to condition the turbine drive fluid in accordance with conditions in 
the real engine. However, the TPCV was a component of the real engine and, as 
such, controlled turbopump speed. 

Chamber pressure (P c ) was computed as a function- of actual 
pump discharge pressure and compared with P c _ re f> which was also generated m 
the computer, and the error, if any existed, actuated the TPCV servo system. 

Since the computed chamber pressure was a direct function of actual pump dis- 
charge pressure which was controlled, indirectly, by the TPCV, the turbine power 
control loop was a direct representation of the corresponding loop in the real 
engine. The Load, Bleed-Line, and Gas-Farm flow-control valves served only 
to transfer computer outputs to the propellant feed system. 

The function of the computer was to relate the load imposed on 
the pump by reactor heating to the energy input at the turbine, which was also a 
function of reactor power output. 
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II. SYSTEM DESIGN 

A. SELECTION AND DESIGN OF PHYSICAL COMPONENTS 

As indicated in paragraph I, A, the C-BEST system was revised to 
form the Generation-I Engine Simulator system in order to carry out the con- 
tractual obligation of Task Item 10 of Contract SNP-I. At the time that the 
C-BEST system was designed, the heat- capacity heat exchanger was designed 
as a reflector segment for heating the turbine drive fluid, with the object of 
minimizing the cost of the system. 

When the C-BEST system was reviewed for adaptation to the Genera- 
tion I Simulator, it was realized that a full-scale reflector would provide more 
realistic engine conditions. However, the cost involved ($120,000) and the time 
to fabricate a full-sized aluminum reflector were overriding considerations. In 
addition, the more sophisticated Gener ation-II system was planned for Phase II 
of the NERVA program, and was being fabricated. Consequently, it was decided 
to continue with the C-BEST system, but to improve its instrumentation system 
and to increase the safety provisions. Some of these modifications have been 
previously listed, and other design considerations are included in the following 
component descriptions 1 

1 . Pump Suction Line 


In the C-BEST system, the existing suction line at the test 
facility was approximately 20 feet long and designed for high pressure, therefore, 
it contained a considerable mass of material. In addition, the tank shutoff valve 
was located about 10 feet from the .pump inlet. These factors caused an unreal- 
istically large heat input to the liquid hydrogen during startup, since the mass of 
material between the shutoff valve and the pump inlet could not be pre-chilled 
prior to startup without also pre-chillmg the pump. 

So that this unrealistic heat input problem could be eliminated 
for the contractual tests, a new line and valve combination was designed, but 
costly facility modifications would have been required if a true mock-up of the 
NERVA line and tank shut-off combination were installed. As a compromise, a 
5 -in. butterfly valve was installed approximately 20 in. above the pump inlet - 
this short section of pipe representing the NERVA suction line - and another 
safety shutoff valve was installed at the tank outlet. The use of two valves m 
the suction supply line allowed the line to be cooled prior to a test, and also 
provided a margin of safety in the event one valve failed during an emergency 
shutdown. 
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A 10-m. Potter mass flowmeter was located between the two 
valves. Although equipped to measure mass flow rate, the flowmeter was used to 
measure volumetric flow only. 

2. Turbopump 


The pump and turbine were mounted on a common shaft, using 
oil lubricated bearings, which were heated electrically to prevent lubricant 
freezing. 


An 11.9-in. -dia aluminum impeller, designated as a "second- 
generation" part, was used in this centrifugal pump. It was radial bladed, had 
inducer vanes at the suction end, and a suction diameter of 7-m. 

A three-stage, pressure compounded, full admission turbine 
was used in the Generation-I system. The total throat area of the first-stage 
nozzle was 1.87 m.^, and the mean blade diameter was 10.98 in. 

3. Heat Capacity Heat Exchanger 


Design of a heat capacity heat exchanger, to simulate the 
engine neutron reflector for the Generation-I Engine Simulator, required estab- 
lishing temperatures and pressures of the fluid at the exit of the heat exchanger 
identical, at any time, with similar properties of the fluid at the exit of the 
engine reflector, that is, the instantaneous heat-exchanger pressure drop and 
rate of fluid heating should always be the same as for the reflector. 

a. Heating the Simulated Reflectors 

Since there was no nuclear heat generation m the simu- 
lated reflector, methods of heating the simulated reflector had to be considered, 
all were discarded because of cost and complexity. In addition, a heat-capacity 
type of heat exchanger provided good startup simulation of the effects of the engine 
reflector. The total heat capacity of the beryllium and graphite reflector, between 
ambient air temperature (520 R) and LH^ temperature (50°R), was calculated to 
be 217,640 Btu. Based upon a typical startup schedule (see Ref. 1), it was found 
that the power generated in the reflector was 4.6 Mw (see Ref. 2) 15 seconds after 
the beginning of a 30-sec turbopump startup, or approximately 34, 000 Btu gener- 
ated by nuclear heating during this period. Since this was only 16% of the energy 
stored as heat capacity, it was postulated that, for the first 15 sec of the startup, 
no simulation of the nuclear heating effect was necessary. 
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b. Cost Factors of Simulation 

The cost of beryllium created another problem of simu- 
lation. Sufficient funds were not available to construct a heat exchanger of 
beryllium in the exact shape of the engine reflector. The first decision, there- 
fore, was to build the heat exchanger as a segment of the engine reflector. 

Because only a percentage of the available fluid was used to drive the turbine, it 
was decided to build a heat exchanger to match only that segment of the reflector 
through which the turbine drive fluid passed. At the time the system was orig- 
inally designed, the cold-bleed engine cycle was the basis for scaling the heat- 
exchanger. In that cycle, turbine weight flow was 13 lb/ sec out of a total of 
83 lb/sec. Thus, a heat exchanger in the ratio of 13/83, or approximately 1/6 
scale, was designed. Since the cost of even this 1/6 scale heat exchanger would 
still have been prohibitive if it were constructed of beryllium, it was decided to 
construct the heat exchanger of aluminum. Although originally chosen largely 
because of availability, subsequent investigation has shown that aluminum pro- 
vides the best available simulation of the thermal properties of beryllium. 

c. Comparison of Possible Heat-Exchanger Materials 

Table II-A, 1 shows a comparison of the thermal proper- 
ties of aluminum, beryllium, and some other materials. The heat capacity of 
the aluminum heat-exchanger is 34,260 Btu between 50 R and 520°R, while the 
heat capacity of T3/83 of the reflector is 34,083 Btu, so the available heat capa- 
cities are nearly the same. 

Additional work to analyze the effect of substituting alumi- 
num for beryllium has been accomplished for the Gener ation-II Engine Simulator. 
These results are also applicable to the Generation I system. From Figure II~A> 1, 
which shows heat capacity as a function of temperature relative to 50 R for both 
aluminum and beryllium, it can be seen that equal amounts of cooling will pro- 
duce almost equal material temperatures, down to approximately 250°R. Since 
the heat exchanger is not cooled to this temperature until about 16 seconds after 
flow initiation, the substitution of aluminum for beryllium is valid for this 
interval. 


The problem remains of insuring that the rate of heat re- 
moval, or the heat-transfer rate, is equivalent for both the aluminum heat exchanger 
and the beryllium reflector. In Figure II-A, 2, which compares the computed 
temperatures of equal sections of beryllium and aluminum cooled at the rates 
indicated by the C-BEST test data, temperatures as a function of time are seen 
to be very nearly the same for both materials. This is partially because the 
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thermal properties, as represented by the inverse thermal diffusivity shown m 
Figure II-A,3, and the thermal conductivity shown in Figure II-A,4, are very 
close above 300°R, but more importantly, because the major thermal resistance 
is the fluid heat transfer coefficient which, for equal geometries and wall tem- 
peratures, is the same for both materials. Thus, equal heat extraction will give 
equal material temperatures and, for these equal wall temperatures, equal flow 
geometries will result in equal cooling rates. 

d. Flow Geometry of Heat Exchanger 

The heat-exchanger pressure drop also depends on flow 
geometry and wall temperature, and, to obtain equal pressure drop as well as 
equal heat transfer, the flow geometry of the heat exchanger should be equal to 
that of the reflector. Since full information concerning the reflector geometry 
was not available when the heat exchanger was designed, calculations were per- 
formed to determine a hole size which would give a pressure drop equal to the 
computed reflector pressure drop. Based upon these calculations, 200 3/16- 
m. -dia. holes, as shown in Figure II-A,5, were calculated to be proper for the 
heat capacity heat exchanger. 

Thermal capacitance was provided by 52 1-in. -thick by 
ll.5-in.-dia. aluminum discs, housed in a 5/8-m. thick stainless steel case with 
flanged and bolted heads. The holes were aligned by two rods piercing all plates, 
and were oriented in eight concentric circles, numbered outward from the center. 

Six CC thermocouples, evenly spaced for the length of the 
heat exchanger, extended into six of the 41 holes in the seventh circle of fuel 
passage holes to record hydrogen temperatures. In addition, six similarly 
located CC thermocouples extending toward the same holes were positioned 
0. 05 to 0. 07-in. short of the coolant holes to record the wall temperature. 

4. Mixing Section 


The mixing section was 5-m-dia. by 27-m. pipe, of Schedule 
40 stainless steel. Flow from the heat-capacity heat exchanger was directed 
axially through the mixing section, and the gas -farm supply was introduced into 
the mixing section at right angles near the upstream end. 

5. Turbine-Power- Control and Shutoff Valves 


The TPCV, designed specifically for the NERVA engine, was a 
3-in. butterfly valve, and was hydraulically operated. The detailed specifications 
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of the electro-hydraulic servo-actuator are given m Appendix B. The actuator had 
been used in the C-BEST system, but the TPCV replaced the 4-m. valve previously 
used, which had been designed for the second stage of the Titan missile. 

A turbine- supply shut-off valve, located downstream from the 
TPCV, consisted of a plug which could be driven shut by an electrically detonated 
explosive charge, and was included as precautionary measure to prevent destruc- 
tion of the turbopump m the event of an overspeed. It could close in 0.006 sec. 
after receipt of an electrical signal, either from the overspeed trip system or 
from the control console. 

6 . Flow-Control Valves and Servo-Actuators 

All three of the flow-control valves and their electrohydraulic 
servo-actuators had also been used in the C-BEST system. Detailed specifica- 
tions of the servo-actuators are given in Appendix B. 

The load valve (LV) was a 5-m. XLR-87 (Titan) butterfly valve, 
modified for hydrogen service. 

The bleed-line valve (BLV) was a modified XLR-87, gas- 
generator, double-poppet-type valve. Flow was allowed through both the oxidizer 
and fuel ports of J;he valve. 

The gas-farm valve (GFV) was a modified XLR-9 1 3-in. butter- 
fly valve. 

The details of the actuator-valve driving linkages are given m 
Appendix B. Since the actuators had been designed for proof pressures of 1000 
psi, it was necessary to operate the cylinders with hydraulic supply pressure of 
800 psig, although the servo-valves were rated for 3000 psi. 
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TABLE II- A, 1 

THERMAL PROPERTIES AT 520°R OF MATERIALS 

Inve rse 
Thermal 

Specific _ Thermal _ ,, 

Density _ , Diffusivity 

Heat p Conductivity . 

C P ^ 3 P 3 K , 

Material Btu/lb°R lb/in. Btu/in. °R Btu/in . - sec°R sec/in. 

Aluminum 

0.215 

0.0978 

0.0210 

0.00205 

10.2 

Be ryllium 

.43 

.0658 

.0283 

.00223 

12.7 

Copper 

.092 

.324 

.0298 

.00517 

5.75 

Paraffin 

.69 

.0 32 

.0221 

.0000034 

6580 

LiH 

.98 

.0284 

.0278 

.000099 

280 

Stainless Steel 

. 108 

.287 

.0310 

.000178 

174 

Mild Steel 

. 100 

.283 

.0283 

.000857 

33 

T itanium 

. 125 

. 16 

.0200 

.000141 

142 

White Pine 

.42 

.0144 

.0060 

.0000167 

360 

Maple 

.42 

.0246 

.0103 

.0000202 

394 
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B. REACTOR SIMULATOR 

1 . Design Considerations 

At the beginning of the C-BEST Program, the portable computer 
available for use as a reactor simulator was equipped with 12 operational ampli- 
fiers, and with the limited funds available, it was possible to build 6 more ampli- 
fiers. But of the 18 amplifiers thus available, 4 would have to serve as error 
detectors to provide the signals necessary to actuate the flow-control valve servos. 
Thus, the problem of designing a reactor simulator for this particular series of 
tests resolved itself to the problem of reducing the reactor equations m such a 
way as to permit a solution with but 14 operational amplifiers, while still re- 
taining the essential dynamic characteristics of the simulated engine components. 
Under ordinary circumstances, an analog computer of sufficient size would be set 
up to solve the complete set of nonlinear equations describing the components to 
be electronically simulated. Reduction of these nonlinear equations included their 
linearizing about the design operating point, taking the Laplace Transform, and 
developing transfer functions by algebraic reduction of the operational equations. 

By using RC networks for the input and feedback impedances of the operational 
amplifiers to simulate transfer functions, maximum utilization of the available 
amplifiers was achieved. 

Reactor simulation must, of course, include the reactor and 
turbine power control loops as well as pertinent functions of the Engine Program- 
mer. In the real engine, the rise in reference chamber temperature (or core 
outlet temperature) is programmed over a period of about 30 seconds to limit the 
approximate rate of change of temperature in the core fuel elements to less than 
200° per second. Considering the fact that the facility tank held only enough 
propellant for a total run time of about 30 seconds, it was necessary to program 
the rise of simulated reactor core outlet temperature to a period less than 20 
seconds, so that a reserve could be maintained for shutdown. Thus, the dynamic 
response characteristics of the prototype reactor control system had to be modi- 
fied to allow the Reactor Simulator to respond adequately to steeper ramp inputs. 
Other simplifications, resulting from the necessity of using linearized system 
equations, were made m the prototype reactor controls. 
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2. Nonlinear Reactor Equations 


The equations describing reactor and reflector characteristics 
are based on the schematic diagram of the cold-bleed cycle nuclear engine, 
Figure U-B-2, 1. The parameter schedule accompanying the figure reflects the 
design point performance data available at the beginning of the C-BEST Program. 
Since there was little change in this parameter schedule at the inception of 
Phase I of the NERVA program, these values were retained for the Generation I 
system. 


A comparison of Figures I-C-l and Il-B-2, 1 shows that, of the 
components comprising the test system, only the turbopump and the TPCV were 
parts of the corresponding nuclear engine. Thus, the Reactor Simulator repre- 
sented the remainder of the nuclear engine, including the pressure drop from the 
pump to the thrust chamber, heat transfer to the fluid in the reflector and reactor 
core, neutron kinetics, and reactor and turbine -power control system character- 
istics. 


The general equations which formed the basis for Reactor Simu- 
lator design, excluding those for the reactor control loop and Engine Programmer 
functions, are listed m the first column of Table XI-B-4, 1. The assumptions and 
simplifications in these equations, which were derived for design point conditions, 
are given below. 


a. Pressure Drop 

The pressure drop in a constant area duct is derived from 
the momentum equation for a compressible fluid, which states that the net force 
acting on a mass of fluid m the line equals mass rate of flow times velocity 
change plus the frictional resistance to flow. Substitution for velocity from the 
equation of continuity yields pressure drop m terms of flow rate squared and 
inlet and outlet densities. The density m the term for frictional head loss is a 
mean value. 


b. Heat Transfer 

According to the design-point data given in Figure 
Il-B-2, 1, the propellant enters the reflector m a gaseous state. Hence, the 
heat-transfer equations for both reflector and reactor core were written for 
gas only. 


The heat generated by radiation in the reflector wall, as 
well as m the reactor core, is proportional to neutron flux level, and the difference 
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between quantity of heat generated and quantity of heat transferred to the fluid by 
convection, represents that quantity stored m the wall. Equations C-l and G-l 
of Table II-B-4, 1 were based on the assumption that the temperature gradient in 
the wall is negligible. The heat transferred to the fluid across the boundary 
results in an increase in the enthalpy of the gas. Film coefficients of heat trans- 
fer are functions of flow rate. 

c. Plenum Capacitance 

In Equations D-l and H-l of Table II-B-4, 1, the time rate 
of change of mass, equal to the net flow, was replaced by substituting, for the 
derivative of gas density, an expression in terms of dP/dt obtained by differen- 
tiating the adiabatic equation, P / = constant. 

d. Turbine Flow Rate 

Because of the lack of computing equipment, it was neces- 
sary to simplify the equation for turbine flow rate by treating the variable area of 
the TPCV as a fixed resistance to flow, which was sufficient for reactor simula- 
tion. 


By equating the equations for flow through the TPCV and the 
turbine nozzle, and assuming ideal throttling occurred, an expression relating 
pressure ratio and area ratio was obtained, 



P 

ro 


a tcv y T cv 
A t K 


( 1 ) 


where 



functions of pressure ratios across the valve and nozzle 


Throughout most of a start transient, sonic flow will exist in the turbine nozzle, 
hence, \}Jt is constant. In the range of valve operation of interest here, the 
pressure ratio across the valve was greater than critical pressure ratio, and a 
plot of area ratio as a function of Pq^/PpQ showed that, m this region, the pres- 
sure ratio changed very little even for large changes of valve area. Thus, 
P Ti/ p r0 ma y be treated as a constant. 
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e. Neutron Kinetics 

The equations for neutron kinetics are well known (Ref- 
erences 3 and 4). Total reactivity was taken as the difference between control 
drum effectiveness and negative temperature reactivity. The positive contribu- 
tion to reactivity by the propellant was not taken into account because it would 
have value only if actual hydrogen density, measured at the heat exchanger exit, 
could be applied to the reactor control loop. Again, the lack of computing 
equipment forced the assumption that, throughout the startup transient, hydrogen 
density at the core entrance was low enough to have negligible effect on reactivity. 
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FULL POWER CONDITIONS: 

OVERALL ENGINE THRUST, LB 56,430 

MAIN THRUST CHAMBER, THRUST, LB 55,000 

TURBINE EXHAUST NOZZLE THRUST, LB 1,430 

OVERALL ENGINE SPECIFIC IMPULSE, SEC 670 

CHAMBER SPECIFIC IMPULSE, SEC 774 

TURBINE NOZZLE SPECIFIC IMPULSE, SEC 110 


PARAMETER SCHEDULE AT FULL POWER 


STATION 

W LB/SEC 

p (pisa) 

T (° R ) 

HEAT, BTU/SEC 
E INPUT 

1 

83 

35 

40 

1 0 

2 

83 

35 

40 

J 

J +40 

3 

83 

800 

50 

4 

83 

700 

120 

J + 35,230 
j + 22,400 

5 

83 

650 

200 

6 

70 

650 

200 

J 0 

J + 1.04 X 10 6 

7 

70 

550 

4090 

8 

70 

550 

4090 

J ° 
) 

9 

13 

650 

200 

J 

10 

13 

500 

200 

J ° 

J -1,950 

11 

13 

50 

140 


COLD-BLEED CYCLE NUCLEAR ENGINE 


SCHEMATIC DIAGRAM 




Fig. II-B-2 ,1 
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3. Specific Reactor Equations 


Values for the system constants appearing m the non-linear 
equations of Table II-B-4, 1 are given in Table II-B-3, 1. The constants and their 
values are divided into two groups one group representing either values taken 
from model specifications for the system shown in Figure II-B-2, 1, or known 
material properties, the second group being quantities derived by substituting 
the values of both known parameters and design point variables into the general 
equations and solving for the unknown constants. Consequently, the substitution 
of the system constants into the general equations yields a set of specific equa- 
tions which will satisfy design point conditions stated m Figure II-B-2, 1. These 
specific equations are shown in the second column of Table II-B-4, 1 and are the 
equations to be linearized. 

The proportionality constants, and (Table III-B-3, 1), 
were derived from a heat balance for design-point conditions, using a value of 
1.22 x 10 15 neutrons/sec for total neutron flux, as based on known core dimen- 
sions and heat flux._ To derive the heat transfer coefficient for the core (h„), an 

o ^ 

arbitrary value of 4500 R was assigned for the core-wall temperature at design 

point conditions. This temperature was reasonable in itself, but also yielded a 
realistic value for h c< 

The heat-transfer coefficient for the reflector (h r ) was calcu- 
lated by Nusselt's equation at design-point conditions, yielding a value of 208°R 
for reflector-wall temperature. In order to obtain a heat balance for the fluid 
in the reflector at the given design point conditions, it was necessary to assign 

a value of 3. 37 to the specific heat (c__) and a value of 3. 82 to c 

r pr pc 

In the specific heat-transfer equations, the factor hA was 
Q Q S 

replaced by Cw ' , where C was found by equating the two factors and substi- 

tuting known design point values. 


TABLE II-B-3, 1 

SYSTEM CONSTANTS FOR REACTOR SIMULATOR EQUATIONS 
A. Given Values of Constants 


Symbol 

Value 

Units 

Symbol 

V alue 

Units 

A 

204 

2 

in. 

A 

60. 2 

2 

in. 

c 



n 



A, 

19. 6 

2 

m. 

A 

25.2 

2 

in. 

d 


r 
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TABLE II - B - 3 , 1 (cont. ) 


Symbol 

V alue 

Units 

Symbol 

V alue 

Units 

A 

t 

2. 64 

2 

in. 

V 

P 

9 

ft. 3 

A 

sc 

90, 300 

2 

in. 

W 

cw 

2385 

lb 

A 

sr 

28,050 

2 

in. 

W 

rw 

2100 

lb 

c 

cw 

0. 54 

Btu/lb°R 

K 

0. 467 


c 

rw 

0. 12 

Btu/lb°R 


0.484 


D 

c 

0. 468 

in. 

A 

0. 00026 


D d 

5. 0 

in. 

A 

0. 00156 


D 

r 

0. 187 

in. 

A 

0. 00191 


g 

32. 2 

, / 2 
ft/ sec 

A 

0. 001995 


k 

1.4 


A 

0. 000607 


k t 

1. 525 x 10" 6 

<5k/°R 

^6 

0. 000173 


r 

2.5 x 10 

sec 

A 

0.0128 

-1 

sec 

L 

c 

52 

in. 

A 

0.0315 

- 1 

sec 

L d 

5 . 66 

ft 

A 

0. 125 

- 1 

sec 

L 

r 

52 

in. 

A 

0. 325 

-1 

sec 

R 

767 

ft. 1° R 

A 

1.55 

-1 

sec 

V 

c 

4. 65 

ft. 3 

^6 

4. 50 

-1 

sec 
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TABLE Il-B-3, i (cont. ) 


B. Derived Values of Constants 


Symbol 

V alue 

Units 

c 

pc 

3. 82 

Btu/lb °R 

c 

pr 

3. 37 

Btu/lb °R 

f 

c 

0. 1325 


f a 

0. 1204 


f 

r 

0. 00187 


h 

c 

4. 88 x 10" 3 

. 2 o 

Btu/sec in. R 

h 

r 

0. 0166 

, 2 
Btu/sec in. R 

K 1 

1. 836 x 10' 11 

Btu/ sec /neutron 

K 2 

8.53 x 10 _1 ° 

Btu/sec/ neutron 

r TCV 

0. 77 
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4. Linearization of Equations 

The total differential of a function of several variables, f 
(x, y, z), is defined as 


Af = Ax + 4^- Ay + Az i 

0 X oy 0 Z 

which states that the differential of a function is equal to the sum of the partial 
differentials obtained by letting the variables change one at a time. Hence, 
each partial derivative is equal to the slope at a point on each curve relating 
the function with one dependent variable. The point in question is the point 
about which the incremental changes in the variables occur. In this case, the 
design point has been chosen as the point for which the linearization is ,to be 
done, and design point variables can be used to evaluate the partial derivatives 
analytically. 


Application of Equation (2) to the specific equations shown in 
Table II-B-4, 1 yields the linear differential equations (in operational form) as 
listed in the third column, Table II-B-4, 1. 

These equations are valid for use only in an operating region 
close to the design point. However, the analog computer will solve these equa- 
tions continuously from zero to design point. Although a steady-state solution is 
possible, the deviation of the solution from design point will be very large be- 
cause of the nature of the linearizations. 

Some simplification and alteration of these linear equations 
(A-3 through J-3) of Table II-B-4, 1 is possible without loss of essential design- 
point time constants. In performing this operation, the object was to obtain a 
set of equations which would have a steady- state solution at the given design point. 

To provide some reasonable basis for simplifying the linear 
equations, steady-state solutions were obtained with the original set of nonlinear 
equations at low power levels. It was found that a change of density with flow 
rate, during a rise in power, had the effect of making pressure drop proportional 
to flow. Hence, the terms expressing incremental changes of density were dis- 
carded, and the proportionality constants in the pressure-flow relations reduced 
by a half. Small terms were discarded. Where the linearization of a product 
of two variables doubled the effectiveness of the incremental variables, the gain 
factors were cut in half. 
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To achieve a heat balance when using design-point values for 
the incremental variables, the coefficients of the heat transfer equations were 
slightly reduced. Fluid temperature at the reflector entrance was assumed to be 
constant since, m the real engine, this temperature can be expected to initially 
drop at a rapid rate to a level not far below its design-point value. 

The resulting set of linear equations were combined to eliminate 
certain variables that are not needed for engine simulation, and were written in 
transfer function form to facilitate construction of a computer circuit utilizing a 
minimum number of amplifiers. (For convenience, the /^-notation was dropped.) 
These equations are listed in the fourth column of Table II-B-4, 1, Equations 
L- 1 through L-8. 


These equations retain the principal dynamic characteristics 
created by flow and heat transfer in the reflector and reactor core, though 
attention is restricted to the region of the design point. Even so, it can be shown 
that the time constants in the nonlinear case do not vary greatly from this one. 
These equations are also satisfied by design-point values of the variables, a con- 
dition required for steady-state operation of the simulated engine. Most im- 
portantly, the effect of an increase of core exit gas temperature, produced by an 
increase in power level, results in an increase in the pressure drop across the 
core . 


The linearization of the neutron kinetics equations is carried 
out m Appendix A, where it is shown that reactor kinetics can be represented, 
over the range of interest here, by a simple transfer function in the form of 
Equation L-7, of Table II-B-4, 1. The net change of reactivity is shown in 
EquationL-8 of that table . 
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o 


REACTOR SIMULATION EQUATIONS 

Equations derived and used for Reactor Simulator design, exclusive of those for the 
Reactor Control Loop and the Engine Programmer 


SPECIFIC EQUATIONS 
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5. Reactor Control Loop 


The linear equations of Table Il-B-4, 1 represent the reflector- 
reactor characteristics required for engine simulation. However, the equations 
for the reactor control loop and the Engine Programmer must be added to com- 
plete the set of equations whose solution is to be instrumented on the analog com- 
puter (Reactor Simulator). 

The NERVA Engine Programmer is a complicated device; it is 
required to act as a sequencing unit, signal generator, and error detector; and it 
is required to perform all of these functions not only for engine control during 
startup, operation and shutdown, but also for safety and malfunction control. 

For purposes of the Generation I Engine Simulator, a single 
integrator was used as a signal generator to provide the reference inputs of 
chamber temperature and pressure to the simulated engine control systems. The 
chief reason for this decision was simply that no more than one operational am- 
plifier could be spared for the Engine Programmer. _ However, since a live 
reactor was not used, the single malfunction and safety system required applied 
only to the turbopump, and control of this system under malfunction conditions 
was vested in a facility automatic overspeed trip circuit. Furthermore, the 
simplification of the simulation and the few tests (four, in all) scheduled to be run 
with the computer permitted only the simplest of startup programs - ramp inputs. 
In addition, since the Reactor Simulator was to start operating at the same time 
propellant flow was initiated, the sequenced startup of the NERVA reactor to the 
initial low power level did not have to be considered 

The ramp inputs of desired temperature and pressure can be 

expressed by 


T 

c 


reference 


Rt, for t < tj 
4090°R for t > t x 


(3) 


P 

c 


reference 


K 3 Rt, for t < tj 
550 psia, for t 


(4) 


The velocity constant (R) is determined by dividing the desired reference value 
by tj, the time required for the ramp to reach the desired value. The constant 
merely relates the T c reference input to P c reference, and allows the use of 
a single integrator for both functions. 
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In the NERVA engine, the reactor-control loop is designed to 
control the rate of change of core exit gas temperature during startup and shut- 
down, and to regulate this temperature at design point. Control of reactor power 
level will be exercised within the temperature control loop. The actual NERVA 
control system will contain shim and regulating control drums, their electro- 
pneumatic actuators, ion- chamber s , and thermocouples for sensing and feeding 
back neutron flux level and gas temperature, respectively. 

Since the representation of such a complex system was not 
warranted for the Reactor Simulator, the simulated control system was designed 
to provide adequate response to a ramp input of T c ~reference. In consideration 
of the requirements that flexibility be provided for selection of the ramp input, 
and that a ramp of 5 seconds might be necessary to make maximum use of the 
limited propellant storage available, a velocity error constant (K v ) of 8. 18 sec"-* 
was taken as the criterion for designing the simulated reactor controller. This 
velocity constant would allow an error of 100°R between T c and T c -reference 
during a 5-sec ramp. Another specification was the requirement for a well- 
damped response, that is, for a damping ratio of about 0.6. 

A block diagram of the reactor control loop is shown m Figure 
II-B-5, 1, where the transfer functions were obtained algebraically through elimi- 
nation of P c and P ro in Equations L-4, 5, and 6 of Table II-B-4, 1 by using the 
relation 


P - P = 1. 43 w (5) 

ro c c 


and then solving the resulting equations for T c and T^ in terms of n, w c , and 
T ro . The controller transfer function (G cr ) is the desired equation for com- 
pleting the Reactor Simulator network. 

Assuming for the time being, that the disturbance functions, 
U] and U 2 > are zero, and letting n Q “ 1.22 x 10^5 j the inner loop of the block 
diagram in Figure II-B-5, 1 can be reduced to 

T 

cw (3.56 X 10 7 ) (S + 0. 1) 

3k~ = (S + 0.0674) (S + 0.308) (S + 262) 


A Bode plot of this transfer function shows that it can be approximated very 
closely by a much simpler expression over the frequency range of interest, 
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T 5 

cw _ 6.55x10 

7k ^ (5.8S t 1) 
cr 


(7) 


A proportional-plus-integral controller will provide the desired response char- 
acteristics. The final control loop configuration is shown in Figure II-B-5.2, 


The controller gain constant (K cr ) can be determined from the 
velocity error constant (K v ) 


Therefore , 


K = 8. 18’ 1 = K (6. 65 x 10 5 ) 
v cr 

K s 1.23 x 10”^, </k/degree error 
cr ° 


( 8 ) 

(9) 


Using this value of K , the open-loop response becomes 


Ef 


= 1.41 f 


r. 


cr , 


cr S (S + 0. 1725) 


( 10 ) 


The closed-loop response becomes 


T c _ 2. 82 (S + 0, 5) 

T 1 “ <S + 2.41) (S + 0. 587) 
c 


where 

the prime notation indicates the reference T . 
r c 

If a root locus plot is prepared of Equation 11 (see Figure II-B-5,3), the con- 
troller time constant ( / f CT ) is determined to be 2 seconds. 

The reactor controller equations can now be written. 



(1.23 x 10" 5 ) 


(2S+ 1) 
S 


( 12 ) 
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and 



(13) 


This controller proved to be stable and provided adequate re- 
sponse characteristics when the simulated reactor control loop was subjected to 
disturbances in flow rate (w c ). See Section II-B-7. 


Page II-B-5, 4 


c 


CONFIDEN TIAL R ET , 









UNCLASSIFIED 



BLOCK DIAGRAM OF TEMPERATURE CONTROL 
LOOP IN REACTOR SIMULATOR 


Figure II-B-5,2 

UNCLASSIFIED 






eGQNF tOCNTtAL-^RP ^ 


Report No. 2223 


6 . Reactor Simulator Circuit 


The transfer functions of Equations L-l through L -6 of Table 
II-B-4, 1, and of Equations 3, 4, 12, and 13, are simulated by the computer 
circuit in Figure II-B- 6 , 1. Scale factors, as indicated in the Figure, were 
chosen so that signals from the pressure and temperature transducers could be 
fed directly into the computer. The maximum output of all pressure transducers 
at Test Stand C -6 is 30 millivolts. With a gam of 700 on the Kmtel-114 trans- 
ducer amplifiers, a full-range signal of 2 1 volts was available. However, for 
convenience a scale factor of 0.01 v/psi was chosen, Scale factors of 0.01 volts 
per psi for T ro , T c , and T cw were also chosen for convenience. The gain of 
the signal from the T mso temperature transducer and associated amplifier was 
0.0305 v/°. This was taken into account by Potentiometers #5 and #6 to match 
the value of 0.025 v/° assigned to T ro . 

Amplifier 5S represents the T c re f erence programmer with 
the limit set at 4090°R. In the initial setup for simulated engine tests, Ampli- 
fier 5S began integrating on "Simulator Run". The remainder of the computer 
circuit began operation at FS-1. 

Potentiometer 6 S at Amplifier 12 represents the coefficient K 3 
m Equation 4. Thus, one input to Amplifier 12 is P c re f erence » 

The four amplifiers. No. 4, 6 , 8 , and 12, are the error de- 
tectors for the feedback control loops in the Engine Simulator, and also include 
the compensation required for stability and performance. The analysis of the 
computer- control system and the determination of compensation is given in 
Section III. 
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7. Simulator Responses 


After fabrication, the analog- computer circuit for reactor 
simulation was checked both statically and dynamically. First, the dynamic 
characteristics of the inner loop of Figure ll-B-5, 1 were determined by applying 
a step change of reactivity (<£k cr ) to Amplifier 1 1 in Figure II-B-6, 1. The 
voltages, P c /100, P rQ / 100, and T ro /40, were fixed at the design-point values 
shown on the circuit diagram. The response of T cw to this step input is shown 
in Figure II-B-7, 1. The time constant of this response agrees very well with 
that of Equation^, and the measured steady-state gain is 6.53 x 1CP °R/c$k. 

Figure II-B-7, 2 shows the transient response of the simulated 
reactor loop of Figure II-B-5,2 to a ramp input of T c ' . As before, P ro /100, 
P c /100, and T ro /40 were fixed. The steady-state velocity lag error between 
T c and T c ' (Figure II-B-7, 2) is 100°R, as specified for the velocity error con- 
stant, Ky., which is 8. 18 sec. 

To determine the stability and performance of the Reactor 
Simulator when P c and P ro are functions of time, i. e. , when the disturbance 
functions and U 2 are not zero, the simulator was subjected to a ramp input 
of PfD Wlt h T c equal to 0. The resulting transients are shown in Figure II-B-7, 3. 

Application of PfQ generated downstream pressures, P ro and P c , 
and the flow rate w r . Since T c ' 1 0, the temperature control loop constrained 
T c to be zero, which means that neutron flux level (n) remained at its initially 
low level (1.8 volts). Consequently, very little change can be detected in T ro . 

About five and one -half seconds after t = 0, a ramp input of 
T c ' was applied while Pf£) was still increasing. The generation of power now 
produced an increase in reflector-gas exit temperature (T ro ). As core-exit gas 
temperature (T c ) and (to a minor extent) T ro increased, chamber pressure (P c ) 
and reflector outlebpressure (P ro ) increased at a more rapid rate so that total 
reflector flow rate (w r ) increased less rapidly than before. Although these effects 
are not quantitatively accurate (because of the linear approximations), neverthe- 
less they are at least qualitatively, the very characteristics required of the 
Reactor Simulator. 

The transient responses shown in Figures H-B-7,2 and 3 were 
very useful during the test program as a means of checking the static and dynamic 
accuracy of the simulator. After the analog computer was shipped to the LRP 
test area at Sacramento, and the circuit reconnected, these transient responses 
were duplicated exactly within the limits of accuracy with which oscillographs 
could be read. Throughout the test program only one amplifier had to be replaced 
(not however, during a a test) and usage of Figures II-B-7, 2 and 3 facilitated the 
detection of the malfunction and subsequent computer checks. 
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III. COMPUTER-CONTROL SYSTEM ANALYSIS 
A. PERFORMANCE CRITERIA 
1. Basis of Design 


A logical and considered approach to the design of a simulated 
engine would be the synthesis of the transfer-controller (load, bleed-line, and 
gas-farm valve servo systems) characteristics that are necessary for integrating 
the physical components and analog computer, but without interposing dynamic 
effects extraneous to the real engine. From these characteristics, dynamic 
response requirements would be established for the design of the transfer con- 
trol systems. 


Considering the fact that the components in the Generation-I 
Simulator, outside of the turbopump, represented a simplified approximation to 
a cold-bleed nuclear rocket engine, that the servo-actuators were already a part 
of the system, and that the analog computer consisted of only 18 operational 
amplifiers, such a synthesis was not applicable. Therefore, the performance 
criteria for the computer-control system, including the reactor simulator circuit 
(design of which is detailed m the preceding section)., were based solely on the 
assumptions that both the hardware and computer were representative of a nuclear 
engine and that, if no steady-state error existed between computer -demanded 
variables and system variables , the dynamic characteristics of this nuclear engine 
would be preserved. The latter assumption implies the use of velocity-type con- 
trollers in the load, bleed-line, and gas-farm valve control loops. It may also 
be inferred from the foregoing assumptions that the turbine-power control system 
should be designed for this representative engine m the same sense that the 
reactor control loop was designed (see II, B, 5). 

Before establishing dynamic response requirements, it is well 
to consider in detail the operation of the transfer controllers. 

2.„ Criteria for Load Valve 


The function of the load valve is to simulate the resistance to 
flow of the reactor core. This is accomplished by comparing computed reflector- 
flow-rate with actual pump delivery, any error actuates the Load-Valve controller 
to move the valve m a direction to reduce the error to zero. This method of 
transferring the computer output to the test hardware was chosen to avoid the 
necessity of computing the flow impedance of the Load Valve and the exhaust duct. 
Actual flow rate is proportional to the A? measured across an orifice near the 
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pump exit. Computed flow rate is a function of pump discharge pressure and was 
converted to a ^P-demand for direct comparison with the measured Ap. Since 
computed flow rate follows the pump pressure, the output of this control system 
must be capable of following a rate of change of Pf£>, another reason for making 
this a velo city-typ e controller. 

3. Criteria for Gas-Farm Control Loop 

There were two reasons for the selection of turbine supply tem- 
perature (T mso ) as the controlled output of the gas-farm control loop, the input 
being computed reflector outlet temperature (T ro ). A primary objective of the 
Generation I Engine Simulator test program was the demonstration of the feasi- 
bility of heat- capacity bootstrapping, therefore the gas-farm valve was required 
to remain closed until sometime after the heat capacity of the heat exchanger had 
diminished. Energy was to be added from the gaseous hydrogen source only after 
the temperature of turbine bleed fluid had dropped below a predetermined value. 
This condition was met by designing the Generation-I system to simulate the cold- 
bleed engine cycle m which the turbine inlet temperature (equal to reflector outlet 
temperature) is 200 R at design operating conditions. Since the temperature of 
the fluid emerging from the heat-capacity heat-exchanger could be expected to 
decay slowly, a considerable time would elapse before T mso would become less 

than T , which follows a uniform rate-of- change of reactor power. Note that, 
r o q q 

from Figure II-B-7, 3 > T ro starts at 120 R and increases to 200 R at design point. 

By permitting the gas-farm controller to be activated only when T mso decreased 

below T ro , the completion of the heat-capacity portion of the starting transient 

was assured. This provision in itself was also a safety feature, because it also 

assured the existence of a load on the pump before the admission of external 

power. 


4. Criteria for Bleed-Line Flow-Control Valve 


In order to relate the turbine energy input to reactor power, it 

was necessary to make both the pressure and the temperature (P and T ) 

at the TPCV equal to simulated reflector outlet pressure and temperature. 
Assuming the fluid to be static in the plenum above the reactor, temperature and 
pressure would be a direct measure of energy and, in the cold-bleed engine, 
plenum conditions would be the same as those at the TPCV (Figure II-B-2, 1). 
Hence, the bleed-line flow-control valve operated to make P mso equal to the 
computed reflector outlet temperature, P rQ . 

5 . Performance Specifications 


Performance specifications for the transfer control systems 
could only be defined in general terms, because the dynamic characteristics of 
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the given servo-actuators were fixed. Previous sections have pointed out the 
necessity of designing the load, gas-farm, and bleed-line control systems as 
velocity controllers capable of following the velocity outputs of the Reactor 
Simulator. Based upon the test results of C-BEST in March 1961, it was pos- 
sible to predict the most probable rate of increase of pump discharge pressure 
that might occur in the Generation-I Simulator. Given such velocities, and con- 
sidering that.the computed command signals would have the same rates, it was 
also possible to establish velocity error coefficients as a guide to the design of 
each of the controllers, provided that permissible velocity-lag errors were first 
selected. The choice of the allowable errors and the selection of velocity error 
coefficients is discussed m the pertinent sections dealing with the analysis of 
individual control loops. In general, servo equalization was designed to provide 
well-damped responses. 

The criteria for designing the TPCV control loop were based 
on the assumption that the computer-hardware-control system was representa- 
tive of a nuclear engine. The transient performance specifications established 
for this engine required that it be capable of following a ramp input of P c refer- 
ence of 110 psi/sec with a velocity-lag error of 25 psi, that P c be held to the 
set-point value of 550 psia with zero error, and that the maximum permissible 
turbine speed overshoot be 10% or less. 

The ramp input was based upon the desire to bring the engine 
to a steady-state operating point in 5 seconds^ obtaining the longest possible run 
time at steady-state. 

The limitation on overshoot arose from considerations for safe 
operation of the turbopump. It was planned to set the over speed trip sensor at a 
level not far above the maximum speed expected in a test. 

A limitation on peak speed was required to avoid a premature 
shutdown during a start transient. 
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B. ANALYTICAL PROCEDURES 
1. Nonlinear Analysis 


Both a nonlinear and a linear analysis were undertaken to deter- 
mine the transient response characteristics of the computer- control system. The 
former involved the application of an analog computer to the solution of the non- 
linear equations describing the characteristics of the physical components. To 
this analog model of the test system, it was planned to connect the portable com- 
puter representing the reactor, and to analyze and synthesize required controller 
characteristics as well as to predict Generation-I system performance. 

As an aid to the formulation of nonlinear system equations, an 
analysis was performed of data obtained m the last of the March 1961 tests on 
theCold-Bleed Engine Simulator, data which was representative of expected re- 
sults from the Generation I Engine Simulator tests. More importantly, the two 
systems were so much alike that the mathematical model of the latter could be 
verified by comparing analog computer results to test data from the former 
system. This was done and the results of that particular study are reported in 
Section III, C, below. 

2. Linear Analysis 


The linear analysis was originally intended as a supplement to 
the nonlinear analysis. As a supplement, a linear system analysis could provide 
great insight into nonlinear system characteristics, and could guide nonlinear 
control- system synthesis by expressing gains and time constants in terms of 
system parameters. 

Toward this end, the linear analysis was undertaken, but un- 
foreseen circumstances forced the abandonment of the analog computer approach, 
and a return to complete reliance on a transfer-function analysis of the system. 
This decision was based on a number of factors. However, the most important 
consideration was that repeated malfunctions of the differential analyzer, coupled 
with a dearth of maintenance personnel, left too vague a period of time in which 
to complete the analysis before the start of an inflexible test schedule. 

The shortcomings of a linear analysis, however, were miti- 
gated by the fact that control of the test system would not be turned over to the 
computer until stable LH 2 appeared m the pump. It was evident, from C-BEST 
heat-capacity starts, that, when stable liquid appeared, turbopump speed was 
over 10,000 rpm and that pump discharge pressure increased rapidly. Similar 
results could be expected in the Generation-I system. 
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Allowing time for visual observation of temperature and pres- 
sure and manual operation of the "Simulator Run" switch, it was reasonable to 
expect discharge pressure as well as speed to be at a high level when the com- 
puter assumed control. Also anticipated was a necessity for limiting maximum 
turbopump speed to somewhere below 20,000 rpm because of uncertainties re- 
garding pump operation at higher speeds. 

Thus, the Generation I Simulator could be expected to operate 
in a restricted range in which the nonlinear ltie s of the system would not be 
severe. The validity of this line of reasoning was demonstrated by the test 
results. 
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C . NONLINEAR SYSTEM ANALYSIS 

This analysis is based on an analog-computer simulation of the Cold- 
Bleed Engine Simulation Test System, with which heat-capacity start tests were 
run in March 1961. Those tests were run without utilizing either the portable 
analog computer for reactor simulation or any of the flow-control valves, the 
system was operated in an open-loop condition. As shown in Figure I-C , 1, the 
hardware consisted of the LH£ turbopump, piping, flow-measuring orifices, 
heat-capacity heat-exchanger, and mixing section. The load valve was closed, 
the bleed-line valve was fully open, the TPCV was fully open, and the gas-farm 
valve was closed. The equations described below were also necessary for the 
linear system analysis. 

1. General Equations - Test System Hardware 


The component equations given in Table III-C,1 were derived 
partly on theoretical grounds and partly by extrapolating empirical coefficients 
from test data. The data chosen for analysis as the most representative was 
obtained in Test No . 1 . 1 -0 1 - ZXF-005 , March 1961, and some pertinent results 
are shown in Figure III-C*1 . All of the variables shown in Figure I-C»2 were 
measured in that test. In addition to yielding empirical coefficients, the data 
analysis also provided a basis for setting up the theoretical equations. 

i a. Flow Equations 

The flow equations were based on the premise that the 
computer would start closed-loop control of the test system only after LH 2 
appeared. at the pump exit. For this condition, the test data indicated that LH 2 
very quickly filled the piping down to the load valve and the heat-exchanger, 
while gaseous conditions prevailed downstream of these sections. Very good 
correlation between flow rates, in various parts of the test system, was ob- 
tained by applying the hydraulic equations for the flow-metering orifices to the 
test data with flow coefficients from the ASME Fluid Meters Report, Reference 5. 
The test data indicated that, for liquid hydrogen, pipe friction losses were 
negligible in comparison to pressure drops across orifices and other large 
discontinuities. Empirical flow coefficients, given in Table III-C,2,were deter- 
mined by applying the theoretical equations to the test data. The bleed-line 
valve and TPCV were treated as orifices, and flow coefficients were extrapolated 
from the Fluid Meters report. 

b. Heat-Capacity Heat-Exchanger 

The constant in the equation for the convective heat- 
transfer coefficient was extrapolated from a heat balance of the transient data 
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based on numerical integration methods. Equation 1 1 of Table III-CjI, which 
was used for determining the temperature of the fluid leaving the heat -exchange r , 
was correlated with experimental data. The best fit with the data was obtained 
by using the average values for L and c v given in Table III-C.2. 

c. Turbopump 

The pump torque and head rise equations were obtained 
from predicted characteristics given in a Technical Memorandum prepared by 
Dept. 8130 of LRP, in February 1961. The turbine torque equation was expressed 
in terms of inlet pressure, temperature, and speed, by combining the equations 
for ideal torque, turbine -nozzle spouting velocity, and efficiency. Turbine 
parameters were obtained from a Technical Memorandum from Dept. 8130, 
December i960. Because no actual pump performance characteristics were 
available, Equations (2), (3), and (4) of Table III - C , 1 were checked against the 
test data. Only fair correlation was obtained. 

2 . Analog-Computer Setup 


Interest had been centered on that part of the start transient 
occurring after LH 2 filled the piping downstream of the pump, therefore no 
attempt had been made to write the two-phase flow equations for the beginning 
of the transient. However, for the analog computer to produce a more realistic 
start transient, the variation of pump discharge density with time was determined 
from test data and programmed on the computer by means of an integrator and 
function generator. 

The function- gene rator output is shown in Figure III-C.l, 
which indicates that, at first, density increased slowly from a very low value, 
then rose rapidly to density of LH 2 • In the pump equations, the density was 
allowed to vary only in the terms. T o be consistent, the inlet temperature 
to the heat-exchanger was also programmed with time according to the test 
results . 


In the turbine equation, the pressure ratio (Px e ^Ti) was taken 
to be constant and equal to 0.09. Sonic flow through the turbine nozzle was 
assumed for all conditions, so the compressible flow coefficient (\ff \) was set 
equal to the sonic value. 

3. Computer Results 


With Pj s constant and equal to 40 psig, corresponding to the 
real test condition, the starting transient shown in Figure III-C,2, was obtained 
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with the analog -computer simulation. The experimentally-determined speed and 
pump discharge pressure are superposed for comparison. 

Table III-Ci3 lists the steady-state values of the variables 
obtained at the peak speed corresponding to liquid density in the pump. These 
steady-state values were later used in the linear analysis for evaluating partial 
derivatives, because they not only are a solution to the nonlinear equations, but 
also agree closely with experimental data. 

From the correlation between computer and experimentsl 
start transients, several conclusions can be drawn that will have an important 
bearing on the linear system analysis. 

a. The startup transient was dominated by the dynamic 
characteristics of the turbopump and by the heat transfer to the fluid. The 
latter can be considered in two distinct stages. First, vaporization of the pro- 
pellant in the suction line and pump produced negligible load on the pump, but 
sufficient flow was generated by tank pressure to accelerate the turbine. The 
higher initial acceleration obtained on the analog model resulted from the 
assumption that sonic conditions at the turbine nozzle prevailed at all times. 
During this first stage, very little heat transfer occurred in the heat- capacity 
heat-exchanger. However, as soon as stable LH 2 appeared in the pump, the 
heat-exchanger dominated energy transfer to the fluid. 

b. Given the proper pump discharge pressure, it was not 
surprising that the steady-state values for flow rates and pressures in the sys- 
tem, obtained at peak speed, should agree closely with the corresponding test 
data, since empirical coefficients had been derived from that data. The important 
factor here is that the predicted pump-head-rise characteristic yielded nearly 
correct values of discharge pressure. 

c. In addition, the predicted pump-torque characteristic 
was verified, by the correlation between the start transients, particularly after 
density in the pump became comparable to a liquid state. This conclusion 
follows from the fact that turbopump response is a function of both M t and M , 

but the turbine torque equation was considered valid, at least for the higher nozzle 
pressure ratios. 


d. The computer results as well as the data analysis indicated 
that steady-state equations could be used for the pre s sure -flow equations, and 
calculations showed that the transport lag in the bleed line was much smaller than 
the time constant of either the heat-exchanger or turbopump. Note that a transport 
lag was actually included, as represented by the bleed-line pressure -drop equations 
coupled with the integration of P mso in. Equation (15) of Table III-C,1. 
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TABLE III-C, I 

NONLINEAR SYSTEM EQUATIONS 
COLD-BLEED ENGINE SIMULATION TEST SYSTEM 
(See Figure III-D,1 for Linearized Equations) 

A. TURBOPUMP DYNAMICS 


T dN M 
I t: = M - M 
at t p 


(1) 


B. TURBINE TORQUE 


M t = < 


CLA.C,r 


F 

i _ 

/P .(k-l)/k 
' te | 

h 


l P txl J 


Vt p - 


vsm ytV 


t ti 


^ i jk 


( 2 ) 


mso 


C. PUMP TORQUE 


w. 


M = C p N 2 +C Q w N-C, 
p pi fD p2 f p3 p fD 


(3) 


D. PUMP DISCHARGE PRESSURE 


w 

+ C . p, n N 2 + C c w, N - C 1 


fD fs p4 fD p5 f p6 p 


fD 


(4) 


where 


P„ = constant 
fs 


E. PUMP DISCHARGE ORIFICE 


W f K fd A f Do V 2g P f D ^ P fd ' P LLV^ 144 


( 5 ) 
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TABLE III-C, 1 (cont.) 


F. BLEED- LINE ORIFICE 


W b K BPO ^PO \J 2& P BPO ^BPO " ^ei^ 144 


(6) 


where 


BP 0 ~ H fD 


G. LOAD- VALVE AND BYPASS ORIFICE 




^ K fyo A fyo + ^V A LV' > 


2g P fyo (P LLV - W^ 44 


( 7 ) 


where 


P fyo ~ P fd 


H. EXHAUST LINE PRESSURE LOSS 


P fDI 14 ‘ T K fDI v i 


( 8 ) 


J. MASS CONSERVATION 


W f \ + 


(9) 


K. BLEED- LINE VALVE 


W h ^ KA bv^ V 2S P fD ^ P LLV ‘ P BP0^ 144 


( 10 ) 
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TABLE III-C, 1 (cont.) 


L. HEAT-CAPACITY HEAT EXCHANGER WALL TEMPERATURE 


dT 


c W — - = - h A (T - T ) 
a a dt s a he 


( 11 ) 


where 


T = 520 R at t = 0 
a 


h = f (w b ) 


0.8 


T = — (T + T ) 
he 2 v Hei Heo ; 

T Hei = 210u ^ t ) " 20 tu(t) + 20 tu(t-8) 


M. HEAT EXCHANGER , CONSERVATION OF ENERGY 


144 


_Hei _ Jfeo j = 7T 8 (l + c T - c T u . ) - 778 
\ P Hei p Heo V He ° v Hei w h 


where 


q. = 


hA (T - T, ) , when T, < T 

s a he ’ he a 


. when T, ST 
7 he — a 


p„ =4.2 
Hei 


144 P 


Heo 


Heo 


R T. 


Heo 


( 12 ) 
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TABLE III-C, 1 (cont.) 


N. HEAT EXCHANGER PRESSURE LOSS 


144 K. 


P - P 
Hei Heo 


1 1 


g A" 


P Heo p Hei b 


P. PRESSURE LOSS, HEAT-EXCHANGER EXIT 


w. 


K. 


'Heo mso 2 p. 


Heo 


(13) 


(l4) 


Q. MIXING SECTION, CONSERVATION OF MASS 


dP 


mso k R 


dt 


lWT T mso v "b 


m 


( w + V - W ) 
g t 


R. MIXING SECTION, CONSERVATION OF ENERGY 


dT 


mso 


dt 


V p 
m mso 


w^ T„ + w T - w, T 
b Heo g g t mso 


where 


T = constant 
g 


p 

mso 


144 P 
mso 

R T 


mso 


S. TURBINE SUPPLY ORIFICE 


(15) 


(16) 


P™. = 0.922 P m 
Ti Ts 


(IT) 
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TABLE III-C, 1 (cont.) 


T. TURBINE POWER CONTROL VALVE, PRESSURE LOSS 

f ^CV^ 


P = P 
Ts mso 


(18) 


U. TURBINE NOZZLE FLOW 


w, 


1 - P 


/ 2g 

OJ 

[ k I 

ik-ll 


[ x _ r (k-l)/k 

1 a - p 4 ) 

/RT 
f mso 

L 


1 - P 4 r 2 / k 



1/2 


P Ti 


where 


A fte 

r P, 
ti 


V. GAS-FARM VALVE FLOW 


w 

g 


= C A 
D gv 



( 20 ) 


where 


P^ = regulated pressure at GFV 
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TABLE III - C , 2 

COEFFICIENTS FOR EQUATIONS OF COLD-BLEED ENGINE 
SIMULATOR TEST SYSTEM 


Symbol 

Value 

Units 

Description 

Source 

a bpo 

1.23 

in. 2 

Orifice area 


A fDo 

11.38 

2 

in. 

Orifice area 


A fY o 

5.11 

2 

in . 

Orifice area 


A 

s 

6200 

2 

in. 

Heat transfer surface 
area 


A 

t 

2.64 

2 

in . 

Throat area, turbine 
nozzle 

* 

c 

a 

0.20 

Btu/lb°R 

Specific heat, aluminum 


c 

V 

c d 

1.86 

0.89 

Btu/lb°R 

Specific heat, hydrogen 

Discharge coefficient, 
empirical 


c i 

4.52 


Coefficient of U/C in 
turbine efficiency 
equation 

* 

C 2 

7.54 


Coefficient of (U/C Q )^ in 
turbine efficiency 
equation 

* 

CPi 

2.526 x 10' 5 


Constant in pump 
torque equation 


C P 2 

7. 19 x 10" 3 


Same 


CP 3 

0.115 


Same 


Cp 4 

3.16x10 3 


Constant in pump head- 
nse equation 



Page III-C-9 




Report No. 2223 


RD P a> , 


TABLE III-C, 2 (cont.) 


Symbol 

V alue 

Units 

Description 

C »5 

1.769 x 10‘ 3 


Same 

Cp 6 

0.0327 


Same 

D b 

3.068 

in . 

Inside diameter, turbine 
fluid bleed line 

°il 

3.068 

in. 

Inside diameter, bypass 
pipe 

D f 

5.047 

in. 

Inside diameter, pump 
discharge line 

g 

32.2 

ft/ sec^ 

Gravitational constant 

I 

0.0733 

ft-lb-sec^ 

Moment - of - ine rtia , 
turbopump rotor 

k 

KA bv 

1.4 


Specific heat ratio 

Effective area of Bleed- 
Line valve , see 
Figure III-D.5 

K 

0.604 


Flow coefficient 

BPO 

K 

0.734 


Flow coefficient 

fD 

K (Dl 

0. 12 


Empirical coefficient 

K fYo 

k lv a lv 

0.517 


Empirical flow coef- 
ficient for APfYo 
orifice and pipe 

Effective area of Load 
valve, see Figure 
III-D, 2 


Source 

** 

** 


* 
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TABLE III-C, 

2 (cont.) 


Symbol 

Value 

Units 

Description 

Source 

K i 

2.24 


Empirical coefficient 


K 2 

0.138 

2. 2 . . 3 . 

sec (in. )(ft ) 

Empirical coefficient 


K 3 

0.922 


Empirical coefficient 


L 

120 

Btu/lb 

Latent heat of 
vaporization 


r 

t 

0.388 

ft 

Turbine mean blade 
radius 

* 

R 

767 

ft/°R 

Gas constant 


V 

m 

0.5 

ft 3 

Mixing section volume 


W 

a 

500 

lb 

Weight of aluminum, 
heat-capacity heat 
exchanger 



sjc 

Technical Memorandum, Dept. 8130, LRP, 30 December 1960. 
Technical Memorandum, Dept. 8130, LRP, 22 February 1961. 
Fluid Meters, Their Theory and Application, ASME, 1959. 
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TABLE III-C, 3 

STEADY-STATE SOLUTION TO NONLINEAR (C-BEST) SYSTEM EQUATIONS 


n t 

= 14,950 rpm 

P Ti 

= 

257 psia 

P £s 

= 54 psia 

P TS 

= 

291 psia 

P fD 

= 407 psia 

T Heo 

= 

358°R 

P 

LLV 

= 386 psia 

T 

TS 

= 

T_ = 368°R 
Ti 

P 

BPO 

= 366 psia 

T 

mso 

= 

365°R 

P 

Hei 

= 348 psia 

T 

Hei 

= 

60°R 

P 

Heo 

= 303 psia 

/> 

mso 

= 

0.156 lb/ft 3 

P 

ms o 

= 303 psia 

W b 

= 

4.21 lb/ sec 


W£ = 48 . 63 lb/ sec 


= 44.2 lb/ sec 
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D. LINEAR SYSTEM ANALYSIS 

In preceding sections, the dynamic characteristics of the Reactor 
Simulator and of test system components were defined. The characteristics of 
the hydraulic flow-control subsystems are analyzed in Appendix B. 

The objective of this section is to determine the kind and degree of 
compensation required in each of the flow-control loops to achieve desired per- 
formance. Neither on the test stand nor on an analog computer is it practical to 
close these loops by trial-and-error methods. The procedure followed here is 
to apply the techniques of linear feedback control system synthesis to determine 
control functions necessary for the achievement of desired performance. 

Of the four feedback control loops comprising the Generation-I 
Engine Simulator, three are necessary for transferring computer outputs to 
the test hardware to simulate nuclear engine performance, and the fourth, the 
turbine -power control loop, represents one of the two real engine control sys- 
tems. The other real engine control loop is incorporated into the reactor 
simulator (portable analog computer). The three subsidiary loops, which control 
the flow of hydrogen through the Load, Bleed- Line, and Gas -Farm valves, were 
analyzed, and the results then incorporated into the analysis of the turbine -power 
control loop; this analysis involves the entire system, including the test hardware 
and the reactor simulator. 

1. Load Valve Control Loop 


The reference input to this control loop is APfp) 0 , generated 
in the reactor simulator as a function of the computed reflector flow rate (w r ) . 
The components in the loop include: the error detecting amplifier (No. 8 in 
Figure II-B-6,1); the Kintel servo amplifier; the hydraulic servo-actuator (with- 
out position feedback); the load valve; the pump discharge line; the APfp> 0 orifice; 
the AP fDo pressure transducer and amplifier; and the signal transmission line 
back to the error detector. 

In linearizing the equations of Table III-Cj 1, was t a ken to 

be constant and equal to 4.2 lb/ft^, and w-^ and PfQ were considered to be input 
disturbances to this control loop. To account for a possible time lag in the 
exhaust duct, the effect of fluid inertia was added to Equation (8), Table III-C,1. 
The linearized equations, in transfer function form, are shown in the block 
diagram of Figure IH-D»1, where the variables are really incrementals . The 
servo-actuator transfer function is described in Appendix B; the gain of the Load 
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valve (K-j^y) is the appropriate slope on the curve of Figure III-D,2; the dynamic 
characteristics of the pipe are defined by K e and ^ e ; and the output APfp) 0 
(in volts) is taken as the output of the transducer amplifier. According to 
Appendix B, K a , the gain of servo-valve amplifier and torque motor, = 25 
milliamps / volt . However, in Figure III-D>1, the torque-motor gain is shown 
to be 2.5 milliamps / volt and the amplifier gain as K a . Splitting the gain in 
this manner was done for expediency, since all the servo amplifiers had been 
modified (subsequently to the experimental frequency responses) to have gains 
variable from 1 to 5 . 

In order to determine K a G c , the input disturbances of Pfj) 
and w^ are, in the first instance, considered to be zero. To define K e and "7T e , 
several operating points were considered, and it was found that, for the 
operating range of interest, K e and 7T e varied very little from the values given 
in Figure XII-D, 1, in which the values were based on design operating conditions. 

From Figure III-B-7,3, the maximum rate of change of w„ 

O * * 

is 13.2 lb/sec^ for a ramp input of P fD = 90 psi/sec. In the computer, w r 
would be 1.32 volts/ sec. The C-BEST test data, Figure III-C,1, indicates 
that higher rates of change of Pf£> (in the order of 140 psi/sec) may be expected. 
Then 


Maximum w r = (140/90)( 1 . 32) = 2.06 volts/sec (14) 


By the linearized relation, 


AP 


fDo 


= 1.372 w r 


(15) 


it is shown that 


d/dt (AP r „ ) — 2.82 volts/sec (16) 

f Do 

If a velocity error coefficient (K y ) is chosen as equal to 2 sec"'*', then the velocity- 
lag error is 1.41 volts, or in terms of flow rate, about 10 lb/sec. A smaller 
velocity error can be expected in the actual system because, according to Figure • 
III-B-7,3, the ramp input of T c re f erence reduces w r to 5 from 13.2 lb/sec. 
Hence, the loop gain will be based on a K y of 2 sec“l . The open-loop gain of the 
block diagram. Figure III-D>1, is 0.332 K a G c ; therefore, K a G c =6.1 for the 
selected K y . 
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The closed-loop roots of the load-valve control system were 
obtained by the root-locus plot shown in Figure III-D,3. Only gain compensation 
was required in this case. The value of K a = 1 was assigned to the servo-valve 
amplifier, and G c =6.1 was assigned to Amplifier 8 in the reactor- simulator 
circuit . 


2. Bleed-Line Control Loop 


The reference input to this loop is P ro (in volts), and the out- 
put is P mso (in psi), or 0.01 P mso (in volts). Included in the loop are the bleed- 
line valve and actuator, the heat-exchanger, the mixing section, the pressure 
transducer and amplifier, and the error detector. 

To simplify the reduction of linearized equations, P LLV and 
Wg were treated as disturbances, and the temperatures in the heat exchanger were 
assumed to be constant; an assumption which was based on the observation that 
heat-exchanger temperatures decayed slowly. In addition, the chief effect of 
the temperature change was an increase in pressure drop, and the assumption of 
constant temperature did not introduce too large an error. 

After combining and reducing the system equations, the block 
diagram of Figure III~D,4 was constructed. The servo-actuator transfer func- 
tion is described in Appendix B. K^y is the slope of the curve in Figure III-D,5. 
The gain K^J), and time constants , and "£ 3 , define the dynamic charac- 

teristics of the bleed line and mixing section. 

The open-loop gain is 0.84 K^G^. The value of K a ^G c ]- ) was 
determined from the root locus plot, Figure III-D, 6 , for a velocity error 
coefficient of 30 sec”^. As shown in Figure III-B-7,3, the maximum rate of 
change of P rQ was 64 psi/ sec, or 6 . 4 volts / sec; therefore, the velocity lag 
error would be only 2.13 psi. The worst error expected was in the range of 
5 psi. To attain an overall K v of 30 sec - -*-, corresponding to a K^G^ of 
35.7; was assigned a value of 5, and G^ (Amplifier 6 ) a value of 7.15. 
Corresponding to the open loop gain of 30, the closed-loop roots of the charac- 
teristic equation were determined as shown on the root-locus plot. 

3. Gas-Farm Valve Control Loop 


The reference input for this control loop is T rQ (in volts), and 
the controlled output is T mso . The block diagram is shown in Figure III-D,7. 
The characteristics of the servo- actuator are the same as for the Load Valve 
(see App. B) . The gain K of the Gas -Farm valve (GFV) is the lesser slope of 
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the curve in Figure III-D, 8 . The gain K 4 is the change in mixing section tempera- 
ture per unit change of valve area, and was obtained by combining the equations 
for the GFY, mixing section, and turbine supply line. When combining these 
equations, the input from the heat exchanger was considered a disturbance, and 
the resistance to flow offered by the TPCV was assumed constant. The gain 
is the gain of the temperature probe and amplifier. The time constant of 
the temperature probe was estimated as 0 . 1 sec. 

The dominant closed-loop roots considered were complex 
conjugate roots for values of open-loop gain sufficiently high enough to yield 
desired velocity lag errors. For that reason, the open-loop gain that would 
provide a closed-loop damping ratio of 0 . 6 was determined first. The closed- 
loop roots are shown on the root-locus plot of Figure III-D,9, and the corres- 
ponding value of root-locus gain is 153,500. Since the system gain is 0.99 
K a gG C g, the required value of K a gG C g is found by 


_ (,0455)(0,1)(153,500) _ 
ag'-'cg - (0 . 028)(5750)( . 99) 


(16) 


Thus, the velocity-error coefficient is 4.35 for the system of Figure III-D,7. 
The expected velocity of T ro was 1.5 volts/sec; therefore, the velocity-lag 
error (Eg V ) would be 0.345 volts. Given an instrumentation calibration of 
0.0305 volts/degree, the velocity-lag error would be 11.3°R. The error 
detector was assigned a gain of one, and the value 4,4 was assigned to the 
servo-amplifier gain, K a . 

4. Turbine-Power Control Loop 


The reference input to this control loop is P^ (in volts). In 
Figure IXI-B-6,1, P^ — 0.1344 T^ at Amplifier 12, where it is compared to the 
controlled variable P c (in volts). P c is generated directly as a result of actual 
pump discharge pressure (Pp-j) fed into the computer as a voltage from the 
transducer amplifier. The error signal (E^p^y) actuates the servo-actuator 
which operates the TPGV. Control of turbine power effectively controls pump 
discharge pressure. Thus, P c is brought into correspondence with P^ . 

In previous sections, the subsidiary control loops were designed 
as though they were non-interacting. That approach allowed the selection of 
error detector and servo -amplifier gains, so that the dynamic characteristics 
of each of the subsidiary control loops are now defined. Each of the subsidiary 
loops can be considered a dynamic element in the turbine power control loop. 
However, the previously neglected disturbance functions, which are really 
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cross -coupling terms, were taken into account in the reduction of the complete 
system block diagram. The reduction was carried through by both algebraic 
and block diagram reduction techniques, and numerical coefficients were re- 
tained throughout the process. In the reduction of the multiple -loop diagram 
many of the transfer function poles were cancelled by zeros which appeared in 
the forward path because poles existed in the feedback paths of subsidiary loops. 

The resulting block diagram is shown in Figure III-D, 10. 

The transfer function (0 o /i) is that of the TPCV actuator, and is developed in 
Appendix B. The gain of the TPCV is the maximum slope of the effective flow 
area curve in Figure III-D, 11. The dynamics of the turbopump and flow loop 
are characterized by a single first-order lag. Since the damping coefficient 
for the turbopump is a function of pump delivery, the effect of the load flow con- 
trol system is to increase the damping on the turbopump. Thus, the turbopump 
time constant of 0.0775 seconds is the result of combining the turbopump and 
load-valve system in the block diagram reduction. Under steady-flow conditions 
at design point, the turbopump time constant was calculated to be 0.100 seconds. 

The roots of the characteristic equation are determined in 
Figure III-D, 12 for a damping ratio £ = 0.62, corresponding to an overshoot of 
10%. The corresponding open-loop gain is K at G ct = 3, and the velocity error 
coefficient is 10.2 sec" . For a ramp input of 110 psi/sec, the velocity-lag 
error between and P c is about 11 psi. In the computer. Amplifier 12 was 
assigned unity gain, and K a ” 3 was assigned to the servo amplifier. 
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IV. TEST PROCEDURES 

A. COMPUTER-CONTROL SYSTEM CHECKOUT 

After the analog computer (Reactor Simulator) was installed in the 
Control Room at Test Stand C-6, the control loops were individually checked for 
continuity, polarity of transducer signals, and polarity of servo-actuator error 
signals. Since the velocity-control loops could not be closed without hydrogen 
flowing in the system, the servo-actuators were checked as positional con- 
trollers, using the sequencing arrangement shown in Figure I-C,l. Taking the 
P m so 1°°P as a typical example, the pressure transducer was removed from the 
test system and a pressure of 650 psi applied by means of a portable hand- 
loader. The resulting signal was brought down to Amplifier 6, Figure II-B-6,1, 
where the reference voltage, P ro , was applied by an initial-condition Pot. on 
the output of Amplifier 1. By raising and lowering P rQ > the bleed-line valve 
was opened and closed, indicating continuity and correct signal polarity. 

In the case of the T mso loop, the temperature resistance bulb was 
removed from the mixing section and dipped into a Dewar of liquid nitrogen, 
providing a fixed reference temperature of -320°F. By increasing or decreas- 
ing the output of Amplifier 3 above or below a voltage corre sponging to -320°F, 
the gas-farm valve was opened or closed. 

No noise problem was encountered with any of the transducer signals 
despite the necessarily high amplification. However, an error in the pressure 
signals in the order of -3.5% was noted during the checkout procedure. This 
loss was attributed to attenuation in the signal transmission lines. 

The static and dynamic performance of the computer alone was 
checked prior to a test run by applying a 5 -sec ramp input of T c re f erence » 
with P.Q 3 — 8 volts, and comparing the transient response of the circuit to that 
of Figure II-B-7,2. 
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B. TEST PROCEDURES 

The procedure followed in all tests was essentially the same: the 

safety valve {Figure I-C, 1) was opened and the pump supply line allowed to fill 
with LH 2 down to the suction valve near the pump; the tank pressure regulator 
was set for tank pressure, as specified for each test, just before countdown; 
countdown began when the supply line had cooled until liquid phase hydrogen (as 
indicated by pressure and temperature probes) existed throughout; on fire- 
switch closure (FS-1), the suction valve opened and LH 2 began flowing into the 
pump (initially at ambient temperature). With the exception of the load valve, 
initial valve positions were as follows: Bleed-line valve open; TPCV open; 

Gas-Farm valve closed. In the first four tests the load valve was closed, but 
in the last four, various starting positions were used. 

In computer -controlled tests, the same procedure was followed, but 
control of the flow-control valves was transferred from the sequence unit to 
the analog computer only when stable LH 2 appeared at the pump exit. This 
point was indicated by a steady pump discharge temperature, T^, of -418°F 
and a smoothly rising pressure, Pfj). 

Shutdown was initiated either when propellant tank level had dropped 
to 15% of capacity, or when a malfunction occurred. On FS-2, the second set 
of relays in the sequence unit was energized. If the computer was in operation, 
the sequencing unit also removed the actuating error signals, leaving only the 
position command signals shown as V 5 to Vg in Figure I-C,l. These signals 
commanded the BLV, TPCV, and GFV to close, thus stopping flow to the 
turbine. Simultaneously, the load valve was opened to increase the load on the 
pump. In the event any of the valves failed to close, a manual switch was to be 
closed immediately to fire the turbine shut-off valve. 
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V. TEST RUNS 


The first four test runs were strictly heat-capacity start tests; that is, the 
computer was not used for engine simulation. However, operation of the computer 
was checked out during two of these runs by feeding the pupnp discharge pressure 
signal to it. The last four runs were computer -controlled tests. Individual changes 
made in the system from run to run are noted below. 

During all the tests, the overspeed trip switch was set for 18,000 rpm. This 
was an arbitrary limit set by the pump group because of a lack of pump data by 
which impeller thrust loads could be predicted. Hence, it was necessary to lower 
the steady-state operating point of the system by reducing the reference set-points 
in the computer. 

A. TEST NO. 1. 1-01-NXS-001 (Figures V-A, 1, 2, 3, 4, 5, 6) 

Tank pressure was set at 22. 5 psig prior to FS-1. The turbine 
started to accelerate at 1800 rpm/sec about 1.0 sec after the fuel suction valve 
opened, and reached a peak speed of 5750 rpm at FS-1 + 39. 2 sec. P reached 
a peak of about 32 psig. Oscillations in P appeared at 16 sec after FS-1 and 
continued through the peak, with a maximum amplitude of about +5 psi. P^g in- 
creased slightly to 25 psig and also exhibited oscillations of 4*1.5 psi. The fre- 
quency of oscillation was about 4 cps on the discharge side and about 1. 5 cps on 
the suction side of the pump. 

The low discharge pressure was attributed to both the low speed apd 
the low fluid density. The low speed occurred as a result of excessive pressure 
loss, 16 psi at the peak point, in the 2 in. ball-check valve downstream of the 
heat exchanger. 
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B. TEST NO. 1.1- l-NXS-002 (Figures V-B, 1, 2, 3, 4, 5, 6) 

Prior to this test, the ball-check valve was replaced by a section of 
two-inch tube . Tank pressure was again set at 22. 5 psig. Because of freezing, 
the fuel suction valve failed to open for a period of 17. 4 sec after FS-1. When it 
did open, the turbopump accelerated at an average rate of 1000 rpm/sec, reach- 
ing a maximum speed of 10, 800 rpm, with a corresponding P equal to 160 
psig. At the same peak, 36 sec after the valve opened, P£g was 18 psig. The 
flow rate at the peak was 31 lb /sec, but the maximum rate was 36. 6 lb /sec about 
6 seconds later. 

Oscillations appeared in P fD 13 seconds after the suction valve 
opened, increasing from zero to a maximum of +25 psi. On the suction side, 
the maximum amplitude was -)-6 psi. This time, the frequency on both sides of 
the pump was the same: about 1.5 cps. 

On the last surge, turbopump speed peaked at 13, 300 rpm with P fD 
being 35 psig. As P^ rose, speed dropped, but recovered to reach the second 
peak at 10, 800 rpm because of che higher flow rate. The drop across the 2-in. 
tube was 8.7 psi at the peak P fD* 
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C. TEST NO. 1. 1-01-NXS-003 (Figures V-C, 1,2, 3, 4, 5, 6) 

For this test, the 2 -in. tube was replaced by a 3 -in. pipe, and a 
2 -in. flapper -type check valve was installed ahead of the heat exchanger . A 
heating coil was wrapped around the hydraulic actuator of the fuel suction valve 
to prevent it from refreezing. 

Tank pressure was set at 62 psig. The turbopump accelerated at 
3200 rpm/sec, reaching apeak speed of 17,000 rpm at 7.7 sec after valve open- 
ing, and following the first dip. At that peak, PfQ was 470 psig, and P^g was 
39. 5 psig. The first peak speed of about 15, 200 rpm occurred 5. 5 sec after 
valve opening when P^ was 75 psig, at the bottom of an oscillation. Speed 
decreased steadily from the second peak to 5500 rpm when FS-2 initiated shut- 
down. Duration was 28.2 sec. Oscillations occurred again at a frequency of 
3 cps, and built up to a peak amplitude of + 40 psi. 

The flapper -type check valve was retained for the remaining tests 
as its operation appeared to absorb only about 7% of the total pressure drop 
across the bleed line. 
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D. TEST NO. 1. 1-01-NXS-004 

This test was a checkout of the turbine over speed shutdown system. 
The overspeed sensing circuit was set to close the FS-2 at a speed of 10, 800 
rpm. Turbine speed was very slowly increased until, at a recorded speed of 
10,780 rpm, the overspeed trip closed FS-2 and fired the turbine shut-off valve. 
The time lapse between FS-2 and the drop in turbine nozzle pressure was deter- 
mined to be 0. 030 sec. 

The computer -controlled tests were to be run with the over speed trip 
set at 18, 000 rpm. Assuming the turbine accelerated from this level at a rate 
of 160, 000 rpm/sec, it would reach a speed of 22, 800 rpm by the time the tur- 
bine shutoff valve closed. Since the maximum safe pump speed is 36,000 rpm., 
there was ample margin for pump protection. 
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E. TEST NO. 1. 1-01-NXS-005 (Figures V-E, 1 through 9) 

This was the first test in which the computer controlled the system. 

Because of the limitation on maximum pump speed, it was necessary 
to lower the P reference set-point in the computer from the design-point value 
of 550 to 495 psia (4. 95 volts). The value of Pjrp which would yield a steady- 
state value of P c equal to 495 in the computer was 680 psia (6. 8 volts). From the 
pump-torque characteristics, the corresponding speed was estimated at 16,500 
rpm. The T c re ference_ set- P 0 i- n t was left at the design value of 4090°R (40. 9 volts) 
Both T c reference an d P c reference inputs started as ramp inputs at the instant 
control was given to the computer, rising from zero to their set-point values in 
5 seconds. The 5-second ramp was chosen to bring the system to a steady-state 
point as quickly as possible. 

Tank pressure was set at 40 psig. The heat-capacity portion of the 
start was similar to occurrences in previous tests. The computer began 
operating on FS-1, accepting P^, and generating w f , P ro , T ro , and P c . Figure 
V-E, 7 shows three of the computer outputs and the T c reference ramp input. 

At 13,4 seconds after FS-1, stable appeared in the pump, and engine con- 

trol was transferred to the computer by the "Simulator Run" switch. At the 
same instant, both T c and Pcreference ram P inputs were started, at which point 
speed was 10, 800 rpm and Pfp) was 210 psig. Almost immediately thereafter 
speed and discharge pressure began to decrease because of the combined effect 
of two factors. As shown in Figure V-E, 8, the load valve opened because W£ 
was apparently less than w r and the TPCV closed because P c was greater than 
P c reference- -As a result of TPCV closure, P mso exceeded the computed 
value of P rQ , so the Bleed-Line Valve closed. The combined effect of closing 
the TPCV, thus cutting turbine power, and of opening the LV, which increased 
pump load, was to cause a sharp decrease in speed. As the ramp input of 
P c reference overtook P c , the TPCV started to reopen, and turbopump speed 
started increasing. The bleed-line valve had opened earlier because the rapid 
decay of P mS o 6 a d allowed P rQ to exceed P m so- The load valve remained open 
until almost the end of the run. In Figure V-E, 8, note that w r reverses polarity 
at about FS-l-f- 16. 4 sec. 

About 18 seconds after FS-1, a loud explosion was heard, so shut- 
down was immediately initiated. However, it turned out that the spark plug on 
the turbine exhaust stack had failed to ignite the hydrogen gas at FS-1; it was 
merely a delayed ignition. It made little difference in run time as tank capacity 
was already low. No damage to the stand resulted from the explosion. 
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A study of the valve error signals, Figure V-E,8, and valve position 
traces indicated that the control loops were operating properly^ although the 
Gas-Farm valve did not open, since the lowest temperature reached in the mixing 
section was about 285°R. In order for the GFV to open, T mso had to go below 
200°R, the steady- state value for T rQ . Although the load-valve actuator responded 
properly to its error signal, a pose-test investigation revealed that the PfDo 
pickup had been inadvertently placed downstream of the orifice where, because 
of pressure recovery, its output signal underwent a sign reversal. Therefore, 
the feedback signal was added to the command signal from w r , and the LV opened. 

In the computer, the sign reversal of w r arose because the ramp 
input of T c re f er ence f° rce d ? r o to become greater than the input Pfp>. 

The noise that appeared on the output of Amplifier 2, did so only 
during a test. No trouble could be found in che amplifier or computer circuit, 
and no noise existed in the instrumentation circuit during pre-run calibrations 
taken at the computer. Therefore, the noise was attributed to vibration at the 
Pfpj pickup. A 0. 005-p fd-capacitor was inserted in the amplifier feedback to 
suppress the noise level in further tests. 
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F. TEST RUN NO. 1 . 1 -0 1 -NXS - 006B (Figure s V -F , 1 through 9) 

The previous test had shown uhat when the ramp input of P c refe rence 
was started at the time of "Sin Run," pump discharge pressure was already 
at a high level. Correspondingly, P c , generated in the computer, was also at 
a high level. The resulting error signal caused the TPCV to close. 

To avoid the possibility of TPCV closure during a heat -capacity 
start and to allow the turbopump to reach steady-state operation, P / 

C X O X C X C 1 X > — c 

was changed from a ramp, to a fixed set-point, which was lowered to 41Z psia, 
corresponding to a Pfj) °f 5Z5 psig. 

Analysis of previous test data showed that a torque balance was not 
achieved with P c reference at 495 psia; the demand for higher discharge pres- 
sures was incompatible with the low turbine inlet temperature for which the 
system had been designed. Calculations were made to determine conditions 
required for a torque balance, and, for P c re f e rence equivalent to 41Z, a 
turbine inlet temperature of -160°F was required. Therefore, the reference 
temperature (T mso ) for operating the Gas-Farm valve was raised to -160°F by 
biasing potentiometer No. 6, Figure II-B-6, 1. Finally, the T c reference ram P 
input was lengthened to Z0 seconds and started at FS-1 to avoid over-riding Pfj). 
the input to the computer from the pump. 

The ^PfDo pressure transducer was relocated to its proper place. 

A diode limiter was placed around Amplifier Z, Figure II-B-6, 1, to keep w r 
from changing sign. 

Tank pressure was again set at 40 psig. The heat-capacity start 
was similar to those of previous runs, Figure V-F, 1. Control was transferred 
to the computer at FS-l-p 1Z. 76 sec, at which time turbopump speed was 
13, 000 rpm and PfD was 310 psig. 

By contrast wi.h the previous test, the Gas Farm Valve, only, 
opened this time, at Z. 3 seconds after "Sim Run. " The valve position traces 
and computer outputs are shown in Figures V-F, 7 and 8, respectively (in the 
latter figure, the Z0 second T c reference ram P input was limited by the recorder 
because the channel gain had been set too high). It is apparent that the GFV 
control loop was stable for almost the entire run, even though T rngo , the con- 
trolled temperature, was oscillatory. The oscillations apparently were caused 
by T^ eo > which was also oscilladng, but with greater amplitude. About two 
seconds before shutdown, the Gas -Farm control loop broke into a diverging 
oscillation. The data indicated that sonic flow existed at the Gas -Farm valve at 
the time the oscillations started. This, plus the fact that Pqj-{ g (GFV supply 
pressure) was also oscillating at the same frequency, leads to *the conclusion 
that the fault lay with the pressure regulators. 

According to the valve error signals (see Figure V-F, 9) the open or 
closed position of each of the other three valves (as shown in Figure V-F -7) was 
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G. TEST RUN NO 1 1-01-NXS-007 

In the previous run, che load valve had remained closed because 
P fD had not reached a level high enough to generate a value for w f m the com- 
puter which was greater than the actual flow. In an effort to obtain more con- 
trol over pump load to achieve steady-state operation, the load valve was 
opened to 23° prior to the start of the test andjto maintain the same flow area 
as before, a 1 5-in -dia orifice replaced che 2 35-in. -dia orifice in the bypass 
section The starting positions of the other valves were unchanged. A torque 
balance based on che data of the preceding test showed that the pump load 
exceeded the available turbine power, therefore, T mSO reference was moved up 
to - 120°F from -160°F. The ramp input for reactor power remained as before. 

When it appeared, from observation of P^, that stable liquid 
hydrogen existed m the pump, control of die system was given to the computer. 
The Load valve immediately closed, turbopump speed exceeded the overspeed 
trip limit, and the system promptly shut down before any damage could occur 

Although, by observation of rapidly increasing P fD ’ there was an 
mdicauon of stable liquid in the pump, the temperature T fD had not yet dropped 
below critical, consequently, the pump was carrying little load The Load valve 
closed because Ap^ Q , higher for two-phase than for liquid flow, was greater 
than the demanded ^P Load-valve closure had che combined effect of unloading 
the pump still more while causing more fluid to go to the turbine. 
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H. TEST NO. 1. 1-01 -NXS- 008 (Figures V-H,l through 15) 

For this test, tank pressure was again set at 40 psig because pump 
suction pressure had dropped to 12 psig during the run of Test -006. 


The P . was further reduced to 345 psia (3.45 volts) to 

c reference 

achieve a steady-state operating point. The pump discharge pressure that 
would produce a P equal to 345 psia, with a T of 4090°R, was about 500 
psia, by the computer equations. Anew torque -balance calculation showed 
that, when T was equal to T^, , or 340 R, a pump discharge pressure of 
460 psig was required. Therefore, the T , set-point was reduced to 

3000°R (30 volts). The correspondingly lower reactor power level produced 
a smaller (simulated) reflector pressure drop, so the required input to the 
computer was 460. 


The steady-state value of T for these conditions was 180 R. How- 

ro 

ever, the output of the amplifier was biased to maintain the set-point temperature 
for the Gas-Farm control loop at -120°F (340°R). It was feared that the GFV 
would open too soon if the reference temperature was raised any higher. 


The T ramp input was lengthened to 30 sec (100°R/sec) 

p t*gic rcncc 

to slow the rate of decrease of w (volts), when the computer assumed control, 
and, in this manner, rapid closure of the Load valve would be avoided. The 
1.5-m. bypass orifice was retained, but the initial Load-valve position was 
increased to 32° for a total flow area of 6 in. 2, to increase the rate of pump 
cooldown during the heat -capacity start. 


Since an oscillation appeared in the Gas-Farm control loop (see 
Figure V - F, 7 ) near the end of the run of Test -006, the gain of the servo- 
amplifier was reduced by a half to provide a greater margin of stability. 


The heat-capacity start was similar to previous runs, and at 15.2 
sec after FS-1, the computer assumed control. Although the Load valve im- 
mediately started to close, pump speed continued to rise until a power satura- 
tion (with the GFV still closed) occurred at 15.8 sec. The LV closed because 
w^ was 48. 6 lb /sec (see Figure V-H, 3) while computed flow (w ) was 40 lb/ sec 
(see Figure V-H, 13), at 15.2 sec. 


Then at 16. 5 sec, the Gas-Farm valve started to open when T , 

(Figure V-H, 10) reached a temperature of 350°R. Since T had already reached 

r o q 

its steady-state level, the reference temperature was supposedly fixed at 340 R. 
This discrepancy is most probably due to some confusion that existed, before 
the run, over the exact calibration of the temperature probe. A second cali- 
bration curve was given which differed by about 10° in this region from a prior 
curve, and the computer was set-up according to the first calibration. 
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As the GFV came open, turbopump speed began co increase again, 
rising to 17,500 rpm, at which point the TPCV began to close, at 18 9 sec 
It may be seen that P c was not recorded in Figure V -H, 13, however, che value 
of P ro . corresponding to a P c of 345 psia, was determined to be 395 psia on the 
computer, prior to the test run It is evident that the TPCV began to close just 
when P c exceeded the P c set-point. (The preceding is true even though T c was 
still below its steady- state value, this is, the ramp input was still increasing 
The reason is that T c had only a small effect on P c , as can be seen from the 
computer circuit) Within two seconds after it started controlling, the Turbine 
Power control system maintained P c at the set-point value, as inferred from the 
fact that P ro was held to 395 psia (3. 95 volts), as nearly as the oscillograph 
could be read (see Figure V-H-13), 

The bleed-line valve had started to close at 17 7 seconds At that 
moment, a P ms o of 420 psia (or rather, 405 psig) was ^he value actually seen 
at the computer, as no attempt had been made to bias the transducer signals 
by 0 147 voUs The corresponding value of P ro was about 380 psi at 17 7 sec. 
During that portion of the run where P ro was held to 395-400 m the computer, 
the controlled variable, P m so> was held co about 415 psig Since the controllers 
were integrating types, the error signal had to be zero Hence, the difference 
of 15-20 psi between P ro and Pj^go was attributed to attenuation of the P m 
signal in the transmission lines (450 feet) and/or in che circuit between the 
transducer amplifier and the computer 

A similar condition existed with respect to Pfp). In the steady-state 
condition, P ro was 395 psia (see Figure V-H, 13), therefore, P c had to be 
345 psia (3 45 volts) to satisfy the computer equations The corresponding pump 
pressure input to the computer had to be 4. 6 volts (460 psia). The actual 
pressure recorded was 495 psig. Figure V -H, 2 Thus, a loss of 0. 35 volts 
occurred in the transmission line, and/or „he circuit between transducer ampli- 
fier and computer 

The ramp input of T c re f e rence> or of power, had a greater effect 
on P r0 , and therefore on w r , than on P c Thus, as the ramp input continued 
to increase, w r decreased, and the load valve gradually closed From about 
20 to 26 seconds, the correspondence between w r and wf was very good, but w r 
was consistently about 2 5 Ib/sec higher than wf This discrepancy could arise 
from the calculation of W£ (see Par. V,A) or from the manner in which 
AP£Dq reference was generated in the computer 

As the load valve slowly closed, Wf was decreased, lightening the 
load on the pump Consequently, the TPCV was gradually closed by its feed- 
back loop, reducing turbine power to hold PfD’ ar *d P c in the computer, steady. 

It is evident that the BLV control loop became unstable about 8 
seconds after the computer began controlling the test The oscillations were 


-e: 
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reflected in P^, and therefore in turbopump speed. About two seconds before 
shutdown, the gas-farm control loop began to oscillate at about 1.7 cps, and the 
bleed-line control loop began to follow the slower oscillation. A divergent 
instability developed throughout the system, resulting m an over speed trip 
shutdown at 28.2 sec. Coincidentally, the shutdown occurred at the same time 
tank capacity reached minimum safe level of 15%. 
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VI ANALYSIS OF TEST RESULTS 

Only the data of Test No. 1 1-01-NXS-008 (see V,H) was selected for 
complete analysis In all previous tests, at least one temperature probe was 
inoperative, in addition, the erroneous location of AP^j-^ differential pressure 
pickup was not discovered and corrected until after the run of Test -005. In 
the eighth test, pressures and temperatures at all flow -metering orifices, and 
on both sides of the pump and heat-exchanger were recorded. 


A. CORRELATION OF FLOW RATES 


The weight flow rates at each orifice (see Figure 1-C,2) were calcu- 
lated by the hydraulic equation, 


w = KA 


2gpA P 


(17) 


If an orifice flowed gas, the expansion factor was taken into account. Thus, 


w 


= KAY M2gpA: 


where 

Y 




as defined in Reference 5 


(18) 


Table VI -A, 1 lists the flow equations for the orifices, the Load valve, and the 
turbine nozzle. In the case of the Load valve, an empirically determined flow 
coefficient (for water) was used to calculate the flow rates of both gas and 
liquid. The compressible flow relation was used for the turbine nozzle 

The calculated flow rates are given in Table VI -A, 2 for a number 
of points taken during the duration of the run. At each point in time, the recorded 
variables were assumed to be at steady- state . The density values were obtained 
from a hydrogen temperature -entropy diagram, and corresponded to the pres- 
sures and temperatures recorded at each orifice. 

The empirical flow coefficient for the Load valve was plotted as a 
function of butterfly -blade angle, which was determined from the valve position 
trace, as shown in Figure V-H, 13 

Only representative points, sufficient for flow-rate correlation, are 
given in Table V-H, 2 Many more points in time were taken to plot the calcu- 
lated flow rates in Figures V-H, 3, 4, 7, and 12. 
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In an effort to determine the degree of accuracy of the equations 
used co determine weight flows, the total weight flows were compared to the 
sum of the component flows In particular, two comparisons are made m 
Table VI -A, Z W£ is compared to (w, +. Wr o +- w^), and w^ is compared to 

(w^-h Both w £ and w t were consistently lower than the respective 

sum of their components. Over the entire run, on the average, W£ and w t 
were 18 5% and 18.6% lower , re spectively, than the sum of their components 
Over the period from 13 6 to 28 seconds, where conditions were steady-state 
or nearly so, W£ and w t were 8 27% and 11 7% lower, respectively, than the 
sum. 


These errors are due to several causes 

1. While data was recorded by a sampling system which recorded 
data at 50 millisecond intervals, still different groups of instrument readings 
were recorded at different time intervals. In some instances, the pressure 
and temperature at a point m the system were recorded at different time 
intervals , 


2 An error was introduced in the summation of component flow 

rates by neglecting the transport lag between the points m the system where the 
component flows were measured 

3. The larger difference between total flows and the component 
summation resulted from the rapid changes m conditions throughout the system 
at the beginning and end of the run. 

4. The error in the flow rate calculated for each orifice is 
probably on the high side, so that the errors were added m the component 
summation 


5. Since no calibration of the orifices was made with hydrogen, 
it was impossible to separate out errors in the calculations made for individual 
orifices However, it is safe to assume that the flow coefficients given by 
Reference 5 are fairly accurate since these coefficients were correlated on the 
basis of Reynolds Number for a large number of fluids. Then the errors 
appearing in the flow-rate correlation were most probably due to errors arising 
m flow rates of the Load valve and turbine nozzle There was an unknown error 
in determination of the blade angle of the Load valve, and another unknown error 
in using a flow coefficient based on water flow tests. The largest error in w^. 
of the turbine nozzle would be due to an error m the use of a single nozzle throat 
area equivalent to the three tandem nozzles that actually existed 
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TABLE VI- A, 1 

FLOW -METERING ORIFICE EQUATIONS 


NOTES 1. Flow coefficients and expansion factors, unless otherwise 
noted, were obtained from the ASME Fluid Meters Report, 
Reference 5 . 


2. Pressure tap locations for all orifices were at 1 D and 1/2 D, 
ahead and after the orifice, respectively. 

3. Reynolds No. was 10^ or greater at each orifice in all 
flow calculations . 


4. Thin-plate, square-edged orifices were used. 


A. AP fT , ORIFICE 

fDo / 2 x 

'V yk (t L / N 2g Ad ap (do 

where 

K 0.734 for £= 0.75 


Y for liquid = 1 


Y for gas = f( — , k,^] 
1 1 

P P - AP 
2 fD fDo 


1 


fD 


d = 3 . 807 in . 
D = 5 .07 in. 
k = 1 .4 


B . AP Bpo ORIFICE 


w = YK 
b 


?Td 


I 2 * ^ 


AP 

BPO ^ BPO 


(2) 
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whe re 


k - 0.604 for^z 0.407 

Y for liquid = 1 

•P_ 


Y for gas = f( — , k,^f 
1 

P P - AP 
2 . BPO BPO 


1 


BPO 


d - 1 . 25 in . 
D= 3.068 in. 
k = 1 .4 


C. A P r ORIFICE 
fyo 


w 


fyo 


= YK 


C 


7T d 


\| 2g Ayo AP fyo 


(3) 


whe re 


K = 0.6176 iox/3- 0 . 49 

Y for liquid = 1 

/ P 2 > 

Y for gas - f( — , k ,/3 

P 9 P. - A P 
2 _ fyo fyo 


1 fyo 

d = 1 . 5 in . 

D = 3.068 in. 
k = 1.4 


D. 


LOAD VALVE 


w 




K 


w 


AP 

AjLV llv 

62.4 


(4) 
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where 


K = empirical flow coefficient obtained from plot of K versus butterfly 

blade angle, and where flow coefficient K was obtained with water 

w 

AP = P - P 

LLV LLV(i) LLV(o) 

~ LLV = fD 


E. TURBINE NOZZLE 


w -= 

1 ( 


1 -/S 


RT 


Ti 


TI 




(5) 



A - 1 . 87 in . 

t 




t 1.87 


A 

R = 767 


7.4 


= 0.504 


k = 1 .4 
F. A P 


gh^s ORIFICE 


w , = YK 

g ^ 2 


T d 


\| 2g ^gh s AP 


gh 2 s 


whe re 


K * 0.724 for^ 0 733 
Y for gas = f (^p~ > k > 
P - AP 1 

P 2 « h 2 S «V 


( 6 ) 


g h 2 s 
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TABLE VI 


CORRELATION OF F 
TEST NO I I — 01 — 


BLEED 

LINE 

FLOW 



Table VI 


2 


LOW RATES 
NXS-008 


GAS 

FARM 

FLOW 

W . 

gh 2 

LB/SEC 


TURBINE FLOW 


W t 

LB/SEC 


LB/SEC 


W, +W . 

b gh 2 

LB/SEC 


0 

0.05 

INV 

INV 

0 

0.20 

INV 

INV 

0 

0.26 

2.53 

0.24 

0 

0.32 

2.82 

0.32 

0 

0.33 

3.44 

0.32 

0 

0 33 

2.63 

0.33 

0 

0.30 

3.72 

Q.42 

0 

0.30 

2.36 

0.28 

0 

0.37 

3.46 

0.35 

0 

0.31 

2.58 

0.31 

0 

0.36 

3.96 

0.45 

0 

0.33 

2.48 

0.35 

0 

0.35 

4.14 

0.45 

0 

0.34 

2.28 

0.28 

0 

0.38 

4.48 

0.52 

0 

0.36 

2.77 

0.43 

0 

0.51 

5.03 

0.60 

0 

0.47 

2.42 

0.28 

0 

0.54 

11.56 

0.65 

0 

0.38 

6.53 

0.31 

0 

0.69 

10.97 

0.56 

0 

0.40 

1.55 

0.32 

0 

1.68 

41.38 

3.14 

0 

4.38 

50.74 

4.96 

2.85 

5.18 

55.65 

5.85 

3.16 

5.29 

49.62 

5.87 

3.28 

5.14 

55.55 

5.46 

2.78 

4.84 

52.93 

5 47 

3.02 

4.42 

52.80 

5.14 

2 62 

3.70 

45.54 

4.36 

2.46 

3.58 

40.87 

4.02 

2.36 

3.76 

44.20 

4.04 

2.93 

3.14 

38.36 

4.05 

1.66 

3.39 

30.83 

3.11 

2.92 

4.69 

56.22 

6.52 

1.52 

3.29 

29.94 

3.97 

2.27 

6.16 

47.98 

7.74 


‘,2 

[FIED 
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B TURBOPUMP PERFORMANCE 


The ideal turbine power calculated by 

w C 

p = - 1 ° 

hdeal ^ 


whe re 


2gk R T 


Ti 


k-1 


(k-lj/k 1 


Te 


Ti 


and is available from Table VI-A.2. 


(19) 


The turbine shaft power was obtained by multiplying the ideal power by the 
turbine efficiency, Figure VI-B.l, and the turbine shaft torque is then given by 


M 

t 


1 > 


( s 


N 


ideal 


rad' 

sec 


, ft-lb 


( 20 ) 


In the steady-state, pump input torque would be equal to turbine 
shaft torque. However, the rotor inertia must be taken into account for the 
transient case , Hence, 

dN 

M r M - J- (21) 

p t dt 

whe re ? 

J = 0.292 ft-lb-sec 

dN/dt is measured on the speed-time plot 


Pump input power is then 


P = N M 

^shaft 


The pump efficiency is 


whe re 


P 

^out 


^P = P P 


/P 

out ^shaft 


144 (P 


f_D_ 

Z 9 . 


P f ) w 

fs f 


fD 


( 22 ) 


(23) 


Page VI-B-1 


CONFIDENTIAL- RD 


Report No . 222 3 


Turbine efficiency was determined, at a given point in time, from 
Figure VI-B.l, for a measured value of tMj / C . The periphe ral velocity was 
U = 0.456 N, ft/sec. 

The results of the calculations are given in Table VI-B,1, in which 
each point in time was assumed to be a steady-state point, and the instantaneous 
acceleration was measured at that point on the transient. The negative values of 
pump-shaft torque indicate that the fluid helped to accelerate the turbopump at the 
start of the transient. During the period of two-phase flow in the pump, most 
of the torque developed by the turbine went into turbopump acceleration. 

At the peak of the heat-capacity start, just before the Gas-Farm 
valve started to open, the shaft power developed by the turbine was 2105 horse- 
power, or almost 47% of design point. 

After the extra energy was added from the gas farm, the ideal 
power available to the turbine was 5027 horsepower at the peak point. Nearly 
half of this power was contributed by the gas flowing yhrough the heat exchanger. 
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TURBOPUMP 




TIME 

SPEED 

SECONDS 

n t 

AFTER FS X 

RPM 


IDEAL 

TURBINE 

TURBINE 

POWER 

EFFICIENCY 

P, 

1 IDEAL 
H.P. 



I. 170 
2.047 
5.119 
6.582 
7.167 
7.606 
7.801 
8.044 
8.191 
8.776 
8.873 
9.068 
9.605 
9.897 

10.043 

10.677 

10.970 

II. 214 
11.506 
11.701 
12.237 
12.725 
13 603 
16.089 
17.991 
19.844 
20.088 
20.771 
22.916 
25.014 
25.891 
26.282 
27.062 
27.208 
27.501 
27.940 
28.233 


30 

909 

4683 

6264 

6813 

7389 

7536 

7926 

7854 

8628 

8484 

8916 

9207 

9744 

9513 

10359 

10194 


960.47 

2877.27 

3455.21 

4062.27 

4166.22 

4058.11 

3816.95 

3816.95 

4162.06 
3812.79 
4157.90 
4041 48 

4153.74 

4153.74 

4361.64 

4245.22 

4848.11 
4232.72 

5022.74 
4112.16 
5326 27 
4315 90 
6178 64 

5633.95 

5742.06 

5712.95 

5733.74 
5671.38 

5791.95 

5821.06 

5862.64 

5688.01 

6033.11 
5937.48 
5363.69 

5563.27 
5939.37 


1.3 

48.2 

90.4 

153.6 

166.8 

143.9 

127.3 

127.3 

186.6 

131.2 

181.2 
157 0 

175.8 

170.8 
21C.5 

188.9 

349.1 

193.0 

396.7 

187.1 

570.0 

217.0 

1867.6 
4048.5 

4973.4 

5027.7 

4920.7 

4533.3 

4317.8 

3650.9 

3583.1 

3542.4 

3328.2 

3480.1 

3929.1 

2965.2 

4555.4 


Table 


VI-B.l 
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C. HEAT CAPACITY HEAT EXCHANGER PERFORMANCE 

The applicability of this analysis to NERVA is questionable on a 
quantitative basis because of the error in heat-capacity heat-exchanger simula- 
tion of the NERVA engine reflector. The flow geometry was not the same as is 
planned for the NERVA reflector, therefore (as demonstrated in Section II. B. 3) 
the similitude of pressure drop and heat transfer for the heat exchanger and 
reflector cannot be established. In addition, since the heat-exchanger segment 
to Cold-Bleed Engine Scale, 13/83 (or 15.7%) of the size of the actual reflector 
but only 5.4 % (Heated-Bleed Engine scale) of the flow was passed through it, 
more heat capacity was available per pound of fluid than is available in the 
real engine and a lower pressure drop occurred. 

Test data for the internal wall and gas temperatures seemed to be 

invalid at the individual stations, when comparing one station against another, 

because calculation of film coefficients indicated contradictory heating and 

cooling between consecutive stations. Averages of the T and T curves 

WHE GHE 

as given in Figures V-H,8 and 9, were therefore used to determine average 
film coefficients as a function of time. From these data, the average heat- 
transfer rates, heat flux, and film coefficients were calculated between FS^ 
and FS^ • 

The equations used for data analysis were as follows 

1. Average Heat Transfer Rate, Q, Btu/sec 

Q = W (ho-hi) (24) 

x x 

where 

ho = enthalpy at outlet, Btu/lb 

hi = enthalpy at inlet, Btu/lb 

W = gas flow rate at time x, lb/sec 
x 

Q = heat transfer rate at time x, Btu/sec 
x 

Subscript x denotes time x 

2 

2. Average Heat Flux, Q / A, , Btu/in. sec 

2 x he 

This quantity is obtained by dividing Equation (24) by the 

total-heat transfer area, A = 6100 m.^ 

he 


was 
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2 o 

3 Average Experimental Film Coefficient, h , Btu/sec in. - R 
g 


Q 


X 


gX \e (T WHE " T GHE ) 


(25) 


whe re 


Q ? 

— — =average heat flux, Btu/in. sec 

Ae 


T = average wall temperature, R 

WHE 

o 

T = average gas temperature, R 

GHE 


2 o 


4. Average Theoretical Film Coefficient, h , Btu/sec in. - R 


, .023 

h = 

g* d -2 


W T 
bx GHE X 


0.8 / 0.40.6, 

c k 
P 


fx 




0.4 


fx 


(26) 


whe re 

, . 2 o 

h = average film coefficient at time, x, Btu/sec in. - R 

g x 

D — 3/ 16-in. -dia hole 


W — gas flow rate at time x, lb/sec 
bx 


^ KA \| 2g A PBPOPBPO inllc l uld 


state 


_ YKA \j 2g Ap BPO P BPO 


in gaseous state 


A = (3/16) x 200— 5.52 in. - flow area 


T = average gas temperature. R 

GHE 

T + T 

T _ GHE WHE f _ o 

fx — _ film temperature, R 
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C = specific heat 
P 

K = thermal conductivity 
p = viscos lty 


'G 0.4 0.6 
_P K_ 

u 0.4 


/fx 


transport properties evaluated for film 
temperature at time x 


Figures VI-C.l and 2 contain plots of the calculated results 
obtained for heat-transfer rates, heat fluxes, and film coefficients. 


The fluid conditions at the entrance to the heat exchanger are 
important because of the concern regarding two-phase flow in the reflector. 
Figure VI-C.3 shows the fluid conditions during Test Run -008 on a hydrogen 
T-S diagram. It can be seen that the fluid conditions at no time lie under the 
dome, therefore, two-phase flow did not exist in the heat exchanger. 
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VII CONCLUSIONS AND RECOMMENDATIONS 
A. HEAT-CAPACITY STARTUP 

The major objective of the Generation I Engine Simulator was 
achieved, that is, the turbopump was successfully "bootstrapped" by utilizing 
the thermal energy initially stored in the metal parts of the propellant feed 
system to heat the turbine drive fluid . In Test No. 1.1-01 -NXS -008 , for 
example, 217 hp was available to the turbine at the point where stable liquid- 
phase hydrogen first appeared in the pump. The speed at that point was 12,500 
rpm, so that sufficient energy was available in the rotor to boost pump pressure. 
As flow increased through the pump, the speed slowed because of the greater 
pump load. However, because of the greater flow to the turbine, the available 
power increased to 1867 hp (at 13.6 sec), and continued to increase as heat 
was extracted from the heat-capacity heat exchanger. Sufficient energy was 
acquired by the turbine fluid in the heat exchanger to boost the turbopump to 
16,200 rpm, which corresponded to the power saturation point in the heat- 
capacity start portion of Test Run -008. The power available to the turbine 
at the power saturation point was 4048 hp, and the turbine shaft horsepower 
was 2105. The corresponding temperature rise of the fluid in the heat exchanger 
was 276°. 


The other heat-capacity start tests (that is, the second and third 
tests m the series) showed that after the power saturation point was reached, 
the transient decay of turbopump speed was much slower than the bootstrap 
phase. The slow decay was the result of the turbine fluid continuing to extract 
nearly enough energy from the heat-capacity heat exchanger to balance the load 
on the pump. By contrast, turbopump speed decayed very rapidly in those 
cases where an overspeed trip caused rapid closure of the turbine shutoff 
valve while the pump was delivering a high flow rate. 

The results of the Generation I Engine Simulator heat-capacity start 
tests can be extended to the NERVA engine. There the bootstrap process would 
be extended to the full power level as a result of reactor core heating. Before 
the heat-capacity of the reflector could be exhausted, reactor power would 
have reached such a level as to greatly heat the fluid in the core. The large 
change in fluid density is essentially an increase in the resistance to flow through 
the core, and consequently, more fluid would be forced through the turbine bleed 
line . 


Further heat-capacity start tests are recommended, but these should 
be conducted with the prototype system (excluding the reactor) to provide more 
realistic flow paths and heat capacities. The heat-capacity startup of a proto- 
type system is expected to be similar to that obtained with the Generation-I 
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Simulator. However, the rate of pump cooldown, and the amount of energy- 
reaching the turbine will be significantly different, sufficiently different to 
justify further tests with a prototype system. 

Provision should also be included in the prototype system to investi- 
gate turbine power control during the heat-capacity portion of engine startup. 

In the Generation I system, the TPCV was held wide open during the start 
transient. This permitted a higher pump delivery than is called for in the 
proposed NERVA startup schedule. In the NERVA engine, thrust-chamber 
pressure is the quantity directly controlled by the turbine power control loop, 
but P is a function of flow from the pump. Therefore, further tests should be 

conducted on the prototype system to determine whether P . should be 

r j l j q retcrcncc 

initially fixed at some percentage of the design point value, or whether it should 

follow a program from the start of propellant flow. It appears from the 

Generation-I test results that scheduling P , from the very beginning 

c reference 

of propellant flow would be the better approach, since pump pressure rose too 
rapidly with this factor fixed at the design point. 
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B. PROPELLANT FEED SYSTEM OPERATION 

Much valuable data and operating experience with a full-scale liquid- 
hydrogen turbopump was gained in these tests. The information obtained will be 
useful in both the design of the propellant feed system, and in the analysis and 
design of the NERVA engine control system. 

Thermal shock caused by flowing liquid hydrogen into an ambient 
temperature pump produced adverse effects neither on the pump nor its opera- 
tion. This results allows the tank shutoff valve to be located upstream of the 
pump, where the location of the valve in this position is especially advantageous, 
because the problem of sealing the pump bearings is much less stringent than 
when the pump is initially immersed in propellant. 

The vaporization of liquid hydrogen in the suction line and pump, 
during the cooldown period, was found to be an inverse function of tank pressure: 
the lower the tank pressure, the longer the time required for chilling the 
suction line and pump. Cooldown took 6 seconds for a tank pressure of 63 psig, 
and as long as 28 seconds with tank pressure set at 23 psig. In the three tests 
where tank pressure was set at 40 psig vaporization in the pump lasted for 
an average period of 1 1 seconds. The shorter cooldown periods are, of course, 
a result of higher initial flow rates produced by higher tank pressures. 

Turbopump acceleration was also a function of tank pressure. More 
flow was available at the turbine for higher tank pressures. Acceleration varied 
directly with tank pressure, from 940 rpm/sec with tank pressure at 23 psig, to 
2960 rpm/sec with tank pressure at 63 psig. The mean was 1100 rpm/sec for 
a tank pressure of 40 psig. 

During the period of 2 -phase flow in the suction line of the pump, 
pressure oscillations appeared on the discharge side of the pump and were 
reflected throughout the system. In most of the heat -capacity starts, the fre- 
quency of oscillation varied from 2 to 3 cps . There was apparently no relation 
of the frequency to either tank pressure or downstream flow impedances. How- 
ever, the amplitude of these oscillations was related to tank pressure and to 
pump speed. Since the higher tank pressures produced correspondingly higher 
accelerations of the turbopump, the amplitudes of the pressure pulses varied 
directly with tank pressure and speed. This fact, plus the fact that the fluid 
density on the suction side of the pump alternated from liquid to vapor, indicated 
that slugs of liquid hydrogen were entering the pump. Apparently, these slugs 
were accelerated by the pump, producing the pressure pulses on the discharge 
side, and, because the speed was greater each time a slug entered the pump, 
the amplitude of each succeeding pressure pulse was greater than the preceding 
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pulse. The oscillations in pressure disappeared when liquid-phase hydrogen 
appeared in the pump. In other words, the suction line and pump were cooled 
down to LH 2 temperature, suggesting that the two-phase flow was largely due 
to heat transfer from the walls of the line and pump. After the appearance of 
stable liquid, pump discharge pressure rose smoothly. No oscillations existed 
during pump discharge pressures greater than the critical pressure, after 
stable liquid-phase hydrogen appeared in the pump. 

The preceding observations on two-phase flow must be qualified 
because of the manner in which fluid density was determined. Only pressure 
and temperature were measured at each side of the pump. The density was 
obtained from a tempe rature -entropy diagram for normal hydrogen (NACA, 

1958) for instantaneous values of pressure and temperature. The response of 
both the temperature probe and pressure pickup, at each location, were 
assumed to be much faster than the system transient. Actually, the temperature 
probe has a relatively slow response which varies with temperature and flow 
rate. Furthermore, a single measurement of pressure and of temperature, at 
any given station, is not sufficient to determine the quality of two-phase fluid, 
nor the stream conditions - stratified, homogenous, or slugging. For future 
tests, attention should be given to the selection of instruments, and the number 
required at a given location, to obtain reliable and accurate measurements of 
two-phase flow of hydrogen. 
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C. COMPUTER-CONTROL SYSTEM OPERATION 

The operation and performance of the analog computer in these tests 
proved the feasibility of computer application to development testing of an un- 
fueled nuclear engine. Considering the degree of simulation achieved with the 
small computer, the most that can be said for the Generation-I Engine Simulator 
is that its behavior was similar to that of a nuclear engine. However, from the 
point of view of computer-control system development, the test program is 
considered highly successful because, in only three attempts, satisfactory 
operation of all four control loops was achieved. Of particular significance in 
this respect is the fact that each of the control loops was only closed with 
hydrogen in the system. 

The operation of the turbine power control loop demonstrated that 

the turbopump can be accurately controlled. During the rise of simulated reactor 

power (Test No. 1 . 1 -0 1 -NXS -008) , the pump flow rate was gradually decreased , 

with a corresponding decrease in pump load. Consequently, the turbopump 

tended to speed up, generating a higher discharge pressure and chamber pressure. 

However, when P„ exceeded P , by a small amount, the TPCV operated 

l c reference 

to reduce speed and discharge pressure. The small error in P existed during 
the period that the flow rate was being reduced. This error was on the order 
of 10 psi, or 2.5%, and is the best measurement than can be made on the 
Sanborn oscillographic record of computer outputs. 

The method of controlling the load flow by comparing the computed 
to the measured flow rate worked very satisfactorily. This method has applica- 
tion to future simulated engine tests because it is considered to be the simplest 
and most direct way of matching pipeline impedance with simulated reactor core 
impedance to flow. 

The Gas-Farm control loop operated satisfactorily, and indicated 
that fairly accurate control of temperature can be achieved. However, the 
Bleed-line control loop tended to ''hunt," because the piston in the BLV actuator 
turned out to be single -ended rather than balanced. The servo-actuator in this 
loop should have had positional feedback to compensate for the unbalanced piston. 

The development of a diverging instability near the end of Test Run 
-008 is attributed principally to the nonlinear characteristics of the system. 

The large change in heat-exchanger exit temperature removed the system from 
the operating range for which the linear system analysis was valid. At the 
point where the computer first took control, the operating conditions corres- 
ponded closely to those on which the linear analysis was based, and the system 
was s table . 
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A linear system analysis is not adequate for the design of a simu- 
lated engine system, particularly if the anticipated range of operation is broad, 
and if more than one control loop is involved. The design of future NERVA 
engine simulators should be based on an extensive nonlinear system analysis, 
which can be conducted with the aid of an analog-computer model of the 
complete system. A linear analysis should be conducted concurrently as a 
guide and as a check on the nonlinear computer analysis. 
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D. HEAT-CAPACITY HEAT -EXCHANGER 

In the case of the heat- capacity heat-exchanger, the quantity and use- 
fulness of the data obtained was somewhat disappointing. The manner in which 
the heat-exchanger had been instrumented did net permit rigorous analysis of 
two-phase heat transfer nor a correlation of heat-transfer coefficients in the 
boiling region. However, as pointed out in the data analysis, Section VI, C, 
fairly good correlation of heat-transfer coefficients based on Nusselt's equation 
was obtained. This correlation was very good where fluid conditions in the 
heat exchanger were gaseous, and rather poor for the boiling region. 

A result of some importance that was obtained in the heat-capacity 
start tests was that, when boiling did occur in the heat exchanger, the pressure 
was above critical and temperature below critical, but pressure oscillations 
did not appear. This indicates that, in the real engine, two-phase flow insta- 
bility will not be a problem in the reflector, provided that reflector exit 
pressure is above critical. The test results indicate that pressure oscillations 
may only appear whenever boiling of liquid hydrogen takes place at pressures 
below critical . 
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VIII. NOMENCLATURE 

Note. See Figure I-C,2 Instrumentation Diagram, for subscripts not shown here 


A , A 
c r 


A. 


A , A 
n t 


A , A 
sc s r 


c , c 
cw rw 


c , c 
pc pr 


D ,D 
c r 

D, 


G , G , G , G 
c cb eg ct 

h , h 
c r 


K 1’ K 2 


K , K , K , K 
a ab ag at 


K 


c r 


Cross-sectional area, reactor core or reflector flow 
passage 

Cross-sectional area, pump discharge duct 

Throat area, thrust nozzle or turbine nozzle 

Heat transfer surface area, core, reflector 

Specific heat, core wall, reflector wall 

Specific heat at constant pressure, core, reflector 

Concentration of ith group of delayed neutron 
precursors 

Diameter of a flow passage , core, reflector 
Diameter of pump discharge duct 
Pipe friction factor 
Gravitational constant 

Gain, error detecting amplifier, as defined in text 
Convective heat transfer coefficient, core, reflector 
Ratio of specific heats 

Power generated per neutron per second, reflector, 
core 

Se rvo-amplifier gain, as defined in text 

Controller gain constant, simulated reactor control 
loop 
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K„ 


K 


V 


L , L 
c r 

L , 


M , M 
t P 

n 


N 

P 

P 

c 

P fD 

P n ,P ro 

P ,P 

ti te 
R 

S 

t 

T 

T 


T , T 
cw wr 

T s T 
ri r o 


T . > T 
ti te 


Negative temperature coefficient of reactivity 
Velocity error coefficient as defined in text 
Mean neutron life time 

Length of a flow passage, core, reflector 

Length of pump discharge duct 

Torque developed, turbine, pump 

Total neutrons in the core 

Turbopump speed 

Pressure, absolute 

Thrust chamber pressure 

Pump discharge pressure 

Reflector inlet or outlet pressure 

Turbine inlet or exhaust pressure 

Gas constant 

Complex Laplace operator 
T ime 

Temperature, absolute 
Thrust chamber temperature 
Wall temperature, core, reflector 
Reflector inlet or outlet temperature 
Turbine inlet or exhaust temperature 
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T 

g 

T 

a 


V ,V 
c 


P 


w 

c 


w 

r 


w 

t 


w 


b 


w 

g 


w„ 


W , W 
cw rw 


0 


l 


A 


i 



/> » P 

’ ri ro 


Vt'Tn 

8k 


Sk 

r 



Temperature of gas -farm supply 

Heat-capacity heat exchanger wall temperature 

Volume, thrust chamber, core inlet plenum 

Weight flow rate, core 

Weight flow rate, reflector 

Weight flow rate, turbine 

Weight flow rate, bleed line 

Weight flow rate, gas farm supply 

Weight flow rate, load valve and bypass orifice 

Total weight, core, reflector 

Fraction of neutrons that are the ith group of delayed neutron 
precursors 

Decay constant, ith group of delayed neutron precursors 

Propellant density, thrust chamber 

Propellant density, pump exit 

Propellant Density, reflector inlet or outlet 

Compressible flow factor, turbine nozzle, thrust nozzle 

Net reactivity, reactor core 

Reactivity contribution of control drums 

Time constant 

Controller time constant, simulated reactor control loop 
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APPENDIX A 

LINEARIZATION OF NEUTRON KINETICS EQUATIONS 


The standard neutron kinetics equations, 


6 

dn _ 8k n {3 n , 

dt " ' X* + A i i 

i=1 


(A-l) 


where 


dC 

dt 




i = l. 


.6 


(A-2) 


can be linearized by developing the expressions for the total differentials of 
n and C : 

l 


af <*> 


8k 


11 * I 

= r 2 ^ + F A(&) " f* ** + L A i c i 


(A-3) 


1=1 


p 


at '“A ' # ^ - [Vi ’] 


(A-4) 


where 


i = l , ... .6 


The parameters Sk Q and n Q are the steady-state values from which the mcrementals 

vary; that is, the variables n and 8k have been replaced by (n + An) and 

(Sk + ASk), where the new variables are the mcrementals. To achieve a steady- 

state m a reactor, the effective multiplication factor must equal unity; hence, 

8k = 8k =0. By performing a Laplace transformation of both equations, the 
o 

second can be solved for AC , 
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3 // 


(A-5) 


s ^ - p? 


3 - 


P i 

s + x. 


A 6^6 

s + X, 


2n 


(A-6) 


6 

Since p = ^ p , the terms in brackets can be combined. The final form of the 
i=l 1 

reactor transfer function, with six groups of delayed neutrons, is 


An 

A^kJ = 


n 

o 





(A-T) 


.* 

Using the values for the p and A and X given m Table II-B-3,1, and ex- 
panding the summation within the parentheses, there is obtained 


An = % 1 (S+0.0128)(S+0.0315)(S+0.125)(S+0.325)(S+4.5)(S+1.55) 

A ^ 6 ’ k ^ S (s +268.2 s 5 +i 6 i 6 s +2132.8 s^+568. 4 s 2 + 36 . 4 i s+0.2708) 


(A-8) 


And, after factoring the denominator. 


An n o 1 

(S+0.0128)(S+0.0315)(S+0.125)(S+0.325)(S+4.5)(S+1.55) 

A(5k) " Jl* S 

( S+ 0 . 0818 ) ( S+l . 41 ) ( S+4 . 4) ( S+ 0 . 008 54 ) ( S+ 0 . 2386 ) ( S+ 262 ) 


or 


An fl| n o (78.2S+1)(31.8S+1)(8S+1)(3.08S+1)(0.222S+1)(0.645S+1) 

A(8'k) k ' l* S ( 12 . 23 S+ 1 ) ( 0 . 71 S+ 1 ) ( 0 . 227 S+ 1 ) ( 117S+1 ) ( 4 . 19 S+I ) ( 0 . 00382 S+ 1 J 

(A-10) 
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A good approximation to the last transfer function is 


An 

A[2kT 


(4.21 x 10 




10 S+l 

S (0.00302 S+l) 


(A-ll) 


A comparison of the two transfer functions is shown m Figure App. A-l. 
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APPENDIX E 

ANALYTICAL AND EXPERIMENTAL DETERMINATION OF CONTROLLER TRANSFER FUNCTIONS 
I. INTRODUCTION 

Each of the four flow-control valves in the Generation-I Engine Simulator 
system is positioned by an electrohydraulic servo -actuator which may be considered 
as a single element m each flow-control loop. In order to conduct both linear 
and non-linear analyses of system stability and performance , it is necessary to 
determine the non-linear equations and, from these, the transfer functions that 
describe the dynamic behavior of each servo-actuator. The non-linear equations 
are required to complete the non-linear analog- computer simulation of the system 
while the transfer functions, derived from the linearized versions of the non- 
linear equations, are needed for linear analyses of each control loop. Thus, the 
objective of this analysis is to derive the controller equations, and to verify 
their validity by comparing theoretically derived frequency response characteristics 
with experimentally determined characteristics. 
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II. DESCRIPTION OF EIECTROHYDRAULIC SERVO ACTUATORS 

The operation of a typical electrohydraulic servo-valve and actuator is shown 
m Fig. App. B-l. The amplifier output produces a differential current m the 
torque-motor coils, resulting m a displacement of the servo-valve spool. Output 
flow from the valve, at a fixed valve pressure drop, is proportional to spool dis- 
placement. Depending upon the polarity of the input signal, the output flow is 
applied to the appropriate end of the actuator piston, to drive the load, i.e., the 
flow-control valve. When used as a positional control system, load position 
(that is, butterfly position) is measured by a potentiometer and fed back for 
comparison with the position command signal. Any error signal is amplified to pro- 
vide a current input to the servo-valve torque motor. Without the position feed- 
back, the servo-actuator will operate as a velocity controller. 

The essential difference between any of the actuator systems lies m the 
linkage used to relate actuator-piston displacement to flow-control valve motion. 

The various linkages used m this program are shown m Fig. App. B-2. Both the 
Load and Gas-Farm valves are butterfly valves rotated by a crank-arm mechanism 
connecting the butterfly shaft to the actuator piston. A rack-and -pinion is used 
to convert actuator translation to rotation of the butterfly shaft m the TPCV. 

In the Bleed-Line valve, a rocker arm produces a 1:1 displacement of both actuator 
and valve pintle. It was necessary to retain the return spring mounted m the 
poppet cavity m the Bleed-Line valve, since the spring is essential for proper 
valve operation. However, the springs normally used m the Load, Gas -Farm, and TPC 
valve actuators were removed to increase their speed of response. Detailed design 
information on each actuator system, as furnished by Dept. 8l4o of LRP, is listed m 
Tables App. B-l, 2, 3; and 4, and servo-actuator characteristics are given m 
Figure App. B-3 to 11. 

A hydraulic power supply system, using hydraulic fluid per MIL-H-5606, pro- 
vided a regulated pressure (P g ) of 800 psig at the supply ports of the servo- 
valves. This system was capable of supplying the maximum flow of 12 gpm for simul- 
taneous operation of all four actuators. 
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III. ANALYTICAL DETERMINATION OF FREQUENCY RESPONSE CHARACTERISTICS 

The analytical development of transfer functions considers each controller 
a positional servomechanism as m Figure App. B-l, since the experimental fre- 
quency responses were obtained for both closed and open loop operation. 

A. DERIVATION OF EQUATIONS 

Except for the actuator- to-valve linkage, the general equations are 
identical for each controller. By considering a load force acting on the actuator, 
the effects of load dynamics can be readily taken into account for different 
linkages. 

The error voltage applied to the servo amplifier is the difference 
between a position command signal, V , and a voltage proportional to the output 
position, 9 0 , of the flow- control valve. 

E = V. - K f 0 q (B-l) 

where 

K f = gam of the position feedback potentiometer, m volts 
per unit displacement 


A differential current, i, is generated m the coils of the servo- 
valve torque motor when a voltage equal to K & E is applied. Because of the coil 
inductance, the current lags the applied voltage. Thus, 


where 


TT ^ + i = K E 
"m dt a 


(B-2) 


= — = torque motor time constant, sec 
m R 


R = coil resistance, ohms 
L = coil inductance, henrys 

Ka = gain of servo amplifier and torque motor, milliamps/volt 
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The no-load flow characteristics of the servo valves are shown m 
Figures App. B-3 and -4. Under no-load conditions, the no-load flow, Qq, through 
the servo valve is given by 


Q = K v/T 
o sv v 


(B-3) 


where 

K = static gain of the servo valve, 
sv 

However, since both servo valves had internal load pressure feedback, the net flow, 
Q, is also a function of the load pressure, as shown m Figures App. B-5 and -6, 
which illustrate the load flow characteristics. If P is supply pressure, and P 
the pressure drop across the actuator piston, then 

« - K v a (b - 4) 

If the piston is stationary and a force, F, is applied to it, a dif- 
ferential pressure, P a , will develop on the piston. To establish this differential 
pressure, the fluid is compressed an amount, AX. A hydraulic spring rate may be 
defined by 


1 “ = A p P a < B -5> 

where 

A = piston area (assume equal areas both sides) m.^ 

P 

= hydraulic spring rate, lb/in. 

The displacement, AX, can be expressed in terms of the bulk ; modulus , (3, of the fluid, 

since AV = A AX, 

P 


A V AP 
_ 2 AV 


( B-6) 


Page 4 
Appendix B 


E NTtA1r~^ RD^ * 


e SONF ID ENTIAL R& - 


Report No. 2223 


and 


since AP = P . Thus, 
a 


AX = 


(v/2)P a 



2 

2 p A 
■ P = 
V 


2 A p 

E 

L 


(B-7 ) 


(B-8) 


where one-half the total actuator volume or stroke, is used because the fluid on 
one side of the actuator is compressed. 

Now, if the piston is moving against an external load, F Q , part of the 
fluid entering the actuator must replace that lost to compression (Ref. 6). From 
Equation (B-6) we have 


dV 

dt 



„ dP 
V i a 

2|3 dt 


(B-9) 


The total flow into the actuator is then the sum of the incompressible and com- 
pressible parts 


dX A 2 dP 

q = A — — + _E — 

H p dt dt 


(B-10) 


where 

3 has been replaced from Equation (B-8) 

X e actuator displacement 

6L 

The force applied to the piston, P & , is in equilibrium (by 
D'Alembert's Principle) with the sum of the inertia force, the dissipative force, 
and the external force on the actuator piston. 
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A P 
P a 




ax 


= M + D — - + 

a ,,2 a dt 

at 


K X + F 
a a o 


( B-ll) 


where 

An actuator spring is included for generality, and 
M & = piston mass 
D = damping factor 

cl 

K = spring rate 

cl 

The external force F q is exerted by the load, which m this case is a flow-control 
valve. For expediency, F q will be derived for each valve m the section dealing 
with the development of individual transfer functions. 
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B. TRANSFER FUNCTIONS 

All of the equations developed in the preceding section, except 
Equation (B-4), are linear. If small displacements of the servo-valve spool are 
considered, it is possible to linearize the flow equation. Thus, 


AQ 



Ai + 
= C 



AP 

a 


(B- 12) 


The partial derivatives are the slopes of the curves in Figures App. B-3, -4, -5, 
and -6, taken at appropriate points. Since the serve valve operates about a null 
position, the incrementals m Equation (B-12) are really total variables. Then 
Equation (B-12) can be written as 

Q = K x i - K 2 P a (B-13) 


The minus sign results from the fact that dQ/dP is negative. 

cl 

The open- loop transfer function for the servo valve and actuator may 
now be derived by combining the Laplace transformations of Equations (B-2), (B-10), 
(B-ll), and (B-13). Thus, 


K K \ 


AM 

Ikm a d ^ 


1 K JD A K 

1 a 

E = 

p a q 3 

2a, pa 

s 2 + 

: 

2a. pa 

^ S+l, 
m J 

Ia + 

i A kT, 

\ p “ ] 

A + A p + K 

i P h 1 


KJC 


*1 

fA 

K_1 

- X + 

S + -A 

a 


A 

P J 


(b-i4) 


In the following sections, F will be related to X m terms of 

O Si 

actuator- valve linkages. 

1. Load Valve 

Linear displacement of the actuator results m rotation of the 
5-m. butterfly valve by means of a slider-crank mechanism. A vector diagram can 
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be drawn relating the external force, F Q , applied by the piston, to a force, F , 
acting at the end of and perpendicularly to the butterfly crank arm. The torque 
applied to the butterfly shaft is F r r, r being the crank-arm radius. Graphical 
resolution of the vector diagram for each actuator position yields the relation 
between F q and T Q , the butterfly torque, shown m Figure App. B-T. The relation 
between actuator displacement and butterfly-blade angle, also determined 
graphically, is given m Figure App. B-8. Comparison of the two plots shows that 
the ratio F q /t o is nearly constant over a considerable portion of the angular 
blade displacement. Then, we may consider that, at least for small servo-valve 
displacements, 


F 

= constant = C (B-15) 

To 1 

o 


Assuming viscous friction m the butterfly shaft bearings and 
considering that hydrodynamic torque on the butterfly, if it exists, is nearly 
linear with blade angle, the equation of motion of the butterfly is 

A s2 + \ s + 1 

From Figure App. B-8, 

0 = K X (B-lj) 

o o a 


0 = T 

o o 


(B-16) 


Combining Equations 


(B-15), (B- 16 ), and (B-17), we obtain 

F = ICK S 2 + DC K S + K C K 
o boo boo boo 


X 

a 


(B-18) 


Note that the factor, C K , has the effect of converting the coefficients of 

o o 

Equation (B-l6) to translational units. Replacing the F q in Equation (B-l4) by 
the F^ in Equation (B-l8) yields 
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X 

a 

E 


(T S+l) 
m 




AM 

Ikm a d \ 

1 K„D A K \ 

Y.Y, " 

S 3 + 

f e + -E^ s 2 + 

2 e + A + -ZS- 

S + 2 e 

_*h 

ip *h / 

Ip p j 

A 

P _ 


(B-19) 


Some simplications were made here by letting 

M = M +ICK (B-20) 

e a b o o 

D = D + DC K ( B-2l) 

e a boo 

K = K +KCK (B-22) 

e a b o o 

Equation (B-19) agrees with Equation (A- 12 of Reference 7> obtained for a dynamic 

2 

pressure feedback servo valve, if K e = 0, A^ is divided out, and K^/A^ is equal to 
the constant of Equation (A- 12). 

2. Gas-Farm Valve 

The open- loop transfer function for the GFV is the same as 
Equation (B-19)> since identical actuators were used on both valves. Only the 
butterfly moment- of- inertia and the damping factor are different. 


3 . Turbine Power Control Valve 

In this case, actuator motion is transferred to the butterfly 
shaft by means of a rack and pinion. Hence, C Q m Equation (B-15) is just the 
reciprocal of the pinion pitch radius, r . By the same token, K q m Equation 
(B-17) is equal to C Q . Then Equations (B-20), (B-2l), and (B-22) become 


M e = M a + h /r p 2 (B - 23) 

D e = D a + (B - 2lt) 
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e a 


K = K + K^/r 


(B-25) 


valve . 


where 


The transfer function of Equation (B-19) also applies to this 


4. Bleed-Line Valve 


The external force F n is given by 


M S + D S + K 
P P P 


X = F 
o o 


X = X 
o a 


Then the equivalent mass, damping factor, and spring rate become 


M = M + M 
e a p 


D = D + D 
e a p 


K = K + K 
e a p 

where 


subscript p refers to the poppet. 


(B-26) 


(B-27) 

(B-28) 

(B-29) 


The moment of inertia of the rocket arm has been neglected. In this case, K e is not 
zero, but is equal to Kp. 
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C. THEORETICAL FREQUENCY RESPONSES 
1. Load Valve Servo 

No flow existed m the flow- control valves during experimental 
frequency response tests; therefore, Kg in Equation (B-19) is zero. Then the 
transfer function becomes 


X 

a 

E 


K K A 
a 1 p 

K D + A 2 
- 2 e P 


S (7r &tl) 
m 


+ — — S + 1 1 

, ,2 £J 

n 
n 


where 


and 



rad/sec 


t 


K n K. M + A 2 D 
2 n e p e 


\/m 

K 

K D + A 2 

V e 

h 

2 e p l 


(B-30) 


(B-31) 


(B-32) 


To evaluate u ) and if , a value for the equivalent damping factor, D g , must be 
estimated. Generally speaking, hydraulic servos have very low damping ratios, so 
a value of D g was chosen to make the damping ratio 0.10. The slope, K^, of the 
load-flow characteristic is obtained from Figure App. B-6, assuming small input 
currents. The hydraulic spring rate, K^, is the slope of the upper portion of the 
curve in Figure App. B-9> the experimentally measured bulk spring rate. 

The torque motor time constant,^* was difficult to determine 
exactly. However, these time constants generally lie between 2 and 7 milliseconds. 
A value for 2r m of h ms was chosen for a first approximation. The parameter values 
are : 
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K = 266,000, (lb/in.) 

= 4.62 x 10"^, (lb sec 2 /m.) 

K 2 = 5.9 x 10 (m.'Ylb sec) 

D = 1.47, (lL sec/in.) 
e 

A = 10.3, (in. 2 ) 

P 

Then 

(jJ 

n 

t 


C o = .75 

O 

= 2400 rad/sec 
= 0.1 


It should be noted that the natural frequency is unaffected by the very small 
damping factor, since is so large. Now the non-dimensionalized open-loop 
transfer function can be written. 


G(S) = 

S (0.004 S + l) 


1 



(2400) 2 


2 ( 0 . 1 ) 

2400 


S + 1 


(B-33) 


2. Gas-Farm Valve Servo 

The only difference between this valve and the Load valve is that 
the moment-of-inertia of the butterfly blade and shaft is much smaller in the GFV. 
Thus, for a damping ratio of 0.1, the natural frequency of the GFV actuator and 
load is 4610 rad/sec. Then the non-dimensionalized open- loop transfer function for 
the GFV servo is 
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C-^S) = 


S (0.004 S + l) 


1 



(46io) 2 


2 ( 0 . 1 ) 
(46ioJ 


S + 1 


(b-34) 


3 . Turbine Power Control Valve Servo 

For this valve also, K g = 0, D g is chosen to yield, a damping ratio 
of 0.1, and 7," is again assumed to be 0.004 seconds. The hydraulic spring rate is 
determined to be 99,000 lb/m. from Figure App. B-10, and the remaining parameters 
from Table App. B-3« 

The natural frequency is 4100 rad/sec and the open- loop transfer 
function is, m non-dimensional form, 


g 2 (s) 


s (o.oo4 s + i) 


(4ioo)‘ 


2 ( 0 . 1 ) 
(Jiioo) 


S + 1 


(B-35) 


4. Bleed- Line Valve Servo 

No spring was contained m the actuator cavity, but it was 

necessary to retain the spring acting on the poppets for proper valve operation; 

therefore, K g = = 320 lb/m., by Equation (B-29)* The actuator spring rate is 

determined, from Figure App. B-ll, to be 119,000 lb/in. The damping coefficient, 

D , for assumed ^ = 0.1, was determined to be 0.J5 lb sec/m. The remaining 
e 

parameters are: 


= 0.004, (sec) 

i/ m 

A = 2.95, (in. 2 ) 

P 

M e = 1.129 x 10 5 , (lb sec 2 /m.) 
= 2.4 x 10"^, (ln.^/lb sec) 
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A comparison of the open and closed-loop transfer functions 
based on Equation (B-19), with K g = 320 lb/in. , and with K g = 0, shows that the 
return spring on the poppet has negligible effect. Hence, the transfer function 
for the Bleed-Line valve is 


G,(S) 


S (0.00L S + l) 


1 



(10,280) 2 


2 ( 0 - 1 ) 

10,280 


s + 1 


(B-36) 
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IV. COMPARISON OF THEORETICAL AND CLOSED- LOOP FREQUENCY RESPONSES 

Both closed- loop and open- loop experimental frequency responses were ob- 
tained for each of the servo-actuator systems. In the open- loop tests, the input 
was the current m the servo-valve coils, connected m series; and the output, 
m one set of tests, was the velocity of the poppet (BLV), or the angular velocity 
of the butterfly blade (TPCV, LV) . The output velocity of the poppet m the 
Bleed-Line valve was obtained with a linear velocity transducer attached directly 
to the poppet. An angular velocity transducer was specially built to measure the 
output blade velocity of the Load and Turbine-Power Control valves. Another set 
of open- loop frequency responses was obtained, but, m these, output position was 
recorded rather than velocity. In all the open- loop tests, the input signal was 
- biased continuously to prevent the actuator piston from drifting off-center; the 
degree of biasing was made proportional to piston drift by a potentiometer bridge 
circuit with the arm of one potentiometer attached to the piston rod. Open- loop 
responses were obtained only with the Bleed-Line, Turbine-Power Control, and Load 
valve actuators, and it was assumed that the data obtained for the Load valve was 
applicable to the Gas-Farm valve. 

The closed- loop frequency responses were obtained with the servo actuators 
set up as positional servomechanisms, as shown in Figure App. B-12. The servo 
valve and actuator are represented by a simplified transfer function which can be 
obtained from Equation (B-30) by dropping the second-order term in the denominator, 
and noting that K^Dg, m the numerator, is negligible by comparison with A^. The 
values of K & and used m the closed- loop tests were derived from Dept. 8130, 

LRP. The value of is the slope on the static flow characteristic taken at a 
differential current of 4 ma. Closed- loop frequency responses were obtained only 
on the Bleed- Line, Turbine-Power-Control, and Gas-Farm valve actuators. The closed- 
loop frequency response of the Load-Valve actuator was assumed to be the same as 
that of the Gas-Farm valve actuator, since the actuators, servo valves, and control 
circuitry were the same m both cases. 
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A. LOAD- VALVE AND GAS-FARM VALVE SERVO ACTUATOR 

The LV and GFV actuators are considered simultaneously because open- 
loop responses were obtained only for the LV actuator, while closed- loop responses 
were obtained only for the GFV actuator. The open- loop experimental responses of 
the Load valve are shown m Figures App. B-13 and -lk. It is evident that there 
is no correlation between the experimental response and the theoretical response, 
as given by Equation (B-33)* Since the theoretical open- loop response for the 
Gas-Farm valve actuator, Equation (B-3^0, is nearly identical to that of the Load- 
Valve actuator up to about 1000 rad/sec, it is clear that no correlation exists 
between the open-loop responses of the GFV actuator. 

In order to determine a reasonable closed- loop response, the data in 
Figure App. B-13 were fitted with the following transfer function: 

= S (0.0455 S + 1) (B-37) 


This transfer function was assumed to apply to the Gas-Farm valve servo 
actuator as well. The open loop response data for blade position output were 
chosen for Fitting (B-37) as that data was considered more reliable than was the 
velocity output data. 

Figure App. B-12 indicates a closed- loop response: 


0 

o 

V 


1 



l/K f 

+ ^ S + 1 


(B- 38 ) 


where 




'K K K 
a 1 o 

A 


(57-3) K 


'h 
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Equation (B- 38 ) applies to both the Load, and Gas -Farm valve servo-actuators. Using 
the gain constants of Figure App. B-12, and with ^ equivalent to 0.0455 sec., the 
closed- loop response becomes 



( 12 . 4 )' 


32 

2 (0.883) 

1274 


s + 1 


(B-39) 


The ciosed-loop response is plotted m Figure App. B-15, where it is compared to 
the experimental closed- loop response obtained with the Gas-Farm valve and 
actuator. Even if a value for TT^ of 0.004 sec is used instead of 0.0455 sec, the 
resulting closed- loop response would be virtually identical to the theoretical 
response shown m Figure App. B-15. 

Since the LV and GFV actuators were to be used as velocity controllers 
m the simulated engine, the extrapolated transfer function, Equation (B-37), was 
used for the dynamic analysis of both the Load valve and Gas-Farm valve control 
loops m the Engine Simulator. 
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B. TURBINE-POWER-CONTROL VALVE SERVO 

The theoretical open- loop frequency response is compared, to experi- 
mental data in Figure App. B-l6, where output (blade) position was recorded. The 
scattering of the data for the smaller input is due to the difficulty of measuring, 
on the oscillograph, amplitude, and phase of the output signal. The trend of the 
data indicates the existence of a quadratic factor with a break-point at about 600 
rather than 4100 rad/sec. The open- loop data obtained with output (blade) angular 
velocity showed the existence of two quadratic factors in the same frequency range. 
However, it was believed that one of these factors was contributed by the network 
used for measuring angular velocity, so that data was discarded. In any case, the 
quadratic factors can be neglected so far as the closed loop response is concerned. 

Experimental and corresponding theoretical closed- loop frequency re- 
sponses are shown m Figure App. B-17- The theoretical closed- loop response is 
characterized by two first-order lags, one at 27-3; and the other at 224 rad/sec. 
The theoretical response plotted in Figure App. B-17 was determined for a K of 
1.95 in. /sec, corresponding to a small- amplitude input. Since the correlation 
was good to at least 32 cps, which was well beyond the bandwidth expected m the 
Turbine-power control loop, the simplified transfer function with a tip °f 0.004 
sec was used in the analysis of the Engine Simulator. 
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C. BLEE D- LINE VALVE SERVO 

The experimental points, for input currents of 4 and 10 ma, peak-to- 
peak, are superposed, m Figure App. B-l8, on the open-loop response of Equation 
( B-36) • There appears to be a quadratic factor present m the range of 70-80 cps. 
However, the gam and phase were difficult to determine from the oscillographs, for 
frequencies greater than 100 rad/sec, because of the large attenuation of the 
output, as well as the effects of backlash m the valve mechanism. 

The closed- loop frequency response, based on Equation (B-36), and 
using the values for K and K„ (amplifier and feedback gams, respectively), as 

a I 

given m Figure App. B-19, is 


0.1431 


V 

1 

r s 2 2(0.738) s ,1 


r ... s 2 _ , 2(0.523) 0 " 


(167.5) 2 \ 


(10, 280) 2 10,280 


(B-40) 


The experimental closed- loop response is also given in Figure App. B-19> where it 
is seen that the correlation is very good to 150 rad/sec. If the closed- loop 
response is obtained from the simplified transfer function, with = 0.004 sec, 
then 



0.1431 

s 2 . 2 (-63) 

(192) 2 192 


s + 1 


(B-4l) 


and it is clear that, for all practical purposes, the simple open- loop transfer 
function, characterized by an integration with first-order lag, is adequate for 
synthesizing the Bleed-Line control loop m the Engine Simulator. 
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TABLE APP. B - 1 


REACTOR-LOAD SIMULATOR 

A. FLOW-CONTROL VALVE 

Nomenclature 

Type 

Total Valve Travel 
Flow Area Vs Blade Position 
Blade Position Vs Actuator 
Displacement 

Force'Torque Ratio Vs Actuator 
Displacement 
Crank Arm Length 
Connecting Rod Length 
Moment of Inertia, Blade, Shaft, 
and Arm 

B. ACTUATOR 

Nomenclature 

Piston Area 
S t r oke 

Piston and Rod Weight 
Equivalent Hydraulic Spring Rate 

C. SERVOVALVE 
Nomenclature 
Type 

Static Flow Gain 
Load-Flow Characteristics 
Supply Pressure 
Coil Resistance 
Rated Current 


VALVE AND SERVO ACTUATOR 


Modified 5 in. TCOV, P/N 1-230650 

Butterfly 

90° ± 1-1/2° 

See Figure III-D.2 

See Figure App . B-8 

See Figure App. B-7 
1.5 in . 

2 .44 in. 

-3 2 

63.7 X10 in. -lb-sec 


TCV Actuator (springs removed), 
P/N 1-218330- 10 

10.3 in. ^ (opening) 1 0. 0 in. ‘'(closing 
2. 12 in. 

4.0 lb 

See Fig. App. B-9 


Moog, Servocontr ols , Inc., 

Model M4369 (or 16-100B) 

Flow control with dynamic pressure 
feedback 

See Figure App. B-4 
See Figure App. B-6 
800 psig 
200 _TL each coil 
12.5 rrulliamperes 
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TABLE APP, B-2 

GAS FARM VALVE AND SERVO ACTUATOR 

A. FLOW-CONTROL VALVE 


Nomenclature 

Type 

Total valve travel 
Flow Area vs Blade Position 
Blade Position vs Actuator 
Displacement 

Force'Torque Ratio vs Actuator 
Displacement 
Crank Arm Length 
Connecting Rod Length 
Moment of Inertia, Blade, Shaft, 
and Arm 

B. ACTUATOR 
Nomenclature 
Piston Area 
Stroke 

Piston and Rod Weight 
Equivalent Hydraulic Spring Rate 

C. SERVOVALVE 
Nomenclature 
Type 

Static Flow Gain 
Load-Flow Characteristics 
Supply Pressure 
Coil Resistance 
Rated Current 


Modified 3 in. TCFV, P/N 
1-224470 
Butte rfly 
90° + 1-1/2° 

See Figure III-D.8 

See Figure App. B-8 

See Figure App. B-7 
1.5 in . 

2.44 in . 

3.7 x 10 3 in. -lb-sec^ 


TCV Actuator (springs removed) 
P/N 1-218330-10 
10.3 in . Z 

10.0 in . (closing) 

2.12 in . 

4.0 lb 

See Figure App. B-9 


Moog Servocontrols , Inc., 

Model M4369 (or 16-100B) 

Flow control with dynamic pressure 
feedback 

See Figure App. B-4 
See Figure App. B-6 
800 psig 
200 -TL each coil 
12.5 milliamperes 
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TABLE APP. B-3 

TURBINE POWER CONTROL VALVE AND SERVO ACTUATOR 
FLOW CONTROL VALVE 


Nomenclature 

Type 

Total Valve Travel 
Flow Area vs Blade Position 
Blade Position vs Actuator 
Displacement 
Pinion Pitch Diameter 
Moment-of-Inertia, Blade, Shaft, 
and Pinion 

Equivalent Mass of Blade, Shaft, 
and Pinion referred to Actuator 
Valve Shaft Resistance Torque 

ACTUATOR 

Nome nclature 
Piston Area 
Stroke 

Piston Weight 
O-Ring Friction 

Equivalent Hydraulic Spring Rate 

SERVO VALVE 

Nomenclature 

Type 

Static Flow Gam 
Load-Flow Characteristics 
Supply Pressure 
Coil Resistance 
Rated Current 


AGC Dwg . 251808-9 

Butte rfly 

90° 

See Figure III-D,11 

See Figure App. B-2 

1.75 in. 

2.111 x 10 3 lb-in.-sec^ 

2.76 x 10 3 lb-sec^/in. 
50 in. -lb (no flow) 


HGV Actuator, Dwg. 1-239790 
4.0 in. 2 (open and close) 

1.90 in. 2 
1 . 14 lb 

65 psig at 1000 psig actuator 
pressure 

See Figure App. B-10 


Moog Servocontrols Inc., Model 
2 106A 

Flow Control 

See Figure App. B-3 

See Figure App. B-5 

800 psig 

200 LL each coil 

12.5 milliamperes 
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TABLE APP. B - 4 

BLEED LINE VALVE AND SERVO ACTUATOR 


A. FLOW CONTROL VALVE 
Nomenclature 

Type 

Total Valve Travel 
Return Spring Rate 
Effective Flow Area vs Valve 
Displacement (both ports) 
Spring Load - Valve Closed 
Weight of Pintle Assembly 
(both pintles) 

Moment of Inertia, Rocker Arm 
Rocker Arm Radius (each side) 

B. ACTUATOR 
Nomenclature 


Piston Area 
Stroke 

Piston Weight 

Equivalent Hydraulic Spring Rate 
C. SERVOVALVE 
Nomenclature 
Type 

Static Flow Gain 
Load-Flow Characteristic 
Supply Pressure 
Coil Resistance 
Rated Current 


Doh.ble Booster, GGV , P/N 
1-223330 
Poppet 
0.6 in . 

320 lb/in. (both ports) 

See Figure III-D,5 
2 64 lb 


-4 2 

7.543 x 10 in. -lb-sec 

1.25 in. 


See Same Dwg. as Flow Control 
Valve 

2 

2.95 in. (open and close) 
Limited to 0.6 in. 

0.203 lb 

See Figure App. B-ll 


Moog Se rvocontr ols , Inc., 

Model 2106A 

Flow Control 

See Figure App. B-3 

See Figure App. B-5 

800 psig 

200 A each coil 

12.5 milliamperes 
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UNCLASSIFIED 



REACTOR-LOAD SIMULATOR VALVE AND 
GAS-FARM VALVE 




MECHANISM FOR CONNECTING ACTUATOR 
DISPLACEMENT TO VALVE MOTION 


Fig. App. B-2 

UNCLASSIFIED 




UNCLASSIFIED 



-16 -12 -8 -4 0 4 8 12 16 
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SERIES CURRENT, i, ma 


GAIN PLOT MO 0 G SERVOVALVE MODEL 2106A 
SUPPLY PRESS-800 PSIG NO LOAD FLOW 


Fig. App. B-3 

UNCLASSIFIED 



UNCLASSIFIED 



»- ' i i i i i i I 

-8 -6 -4 -2 0 2 4 6 8 

SERIES CURRENT, i, ma 


GAIN PLOT MOO G SERVOVALVE MODEL M4369 
SUPPLY PRESS-800 PSIG,NO LOAD FLOW 


Fig. App. B-4 
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CODE: 


14 ma 
12 ma 
10 ma 
8 ma 
6 ma 
4 ma 
2 ma 
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LOAD FLOW CHARACTERISTICS 
M00G SERVOVALVE MODEL 2106A 
SUPPLY PRESS -800 PSIG 


Fig. App. B-5 
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CODE: 


12,14,16 ma 
10 ma 
8 ma 
6 ma 
4 ma 
2 ma 
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CL 



0 100 200 300 400 500 600 700 800 
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LOAD FLOW CHARACTERISTICS 
MOOG SERVOVALVE MODEL M4369 
SUPPLY PRESS -800 PSIG 


Fig. App. B-6 
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0 .25 .50 .75 1.00 1.25 1.50 1.75 2.00 2.25 

PISTON DISPLACEMENT, X , INCHES 

a ' 


LOAD AND GAS-FARM VALVES TORQUE RATIO 


Fig. Add. B-7 
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RADIANS 


UNCLASSIFIED 


NOTE: 


6 = K X 
o o 













PISTON DISPLACEMENT, X, , INCHES 

cl 


BLADE POSITION AS A FUNCTION OF ACTUATOR DISPLACEMENT 
FOR LOAD VALVE AND GAS FARM VALVES 


Fig. App. B-8 
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CODE: 

PISTON CENTERED (VALVE 50% OPEN) 

VOLUME @ 0 psig = 948 ml 

PISTON EXTENDED (VALVE OPEN) 

VOLU ME @ 0 psig = 480 ml 

PISTON RETRACTED (VALVE CLOSED) 

VOLUME @ 0 psig = 928 ml 



CHANGE IN VOLUME, AV, MILLILITRES 


BULK SPRING RATE 
GFV, RLSV AND ALT, TPCV 
ACTUATORS 


Fig. App. B-9 
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NOTE- 


PISTON EXTENDED (VALVE OPEN) 

VOLUME® 0 psig = 130 ml 

PISTON RETRACTED (VALVE CLOSED) 

VOLU ME @ 0 psig = 278 ml 

PISTON CENTERED (VALVE 50% OPEN) 

VOLUME® 0 psig = 289 ml 



BULK SPRING RATE PROTOTYPE TPCV 


Fig. App. B-10 
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CODE: 

PISTON CENTERED, VALVE 50% OPEN 

(VOLUME @0 psig = 117 ml) 

PISTON EXTENDED, VALVE OPEN 

(VOLUME @0 psig = 55 ml ) 

PISTON RETRACTED, VALVE CLOSED 

( VOLUME @0 psig = 90 ml) 



0 1.0 2.0 3.0 

CHANGE IN VOLUME, AV, MILLILITRES 


BULK SPRING RATE BLV ACTUATOR 


Fig. App. B-ll 
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LEGEND: 

K = 25 MA/VOLT 

cl 

K x = 2.15 IN. 3 /SEC/MA 
K =0.75 RAD/IN. 



LOAD OR GAS-FARM VALVE SET-UP 
FOR CLOSED-LOOP FREQUENCY RESPONSE 


Fig. App. B-12 
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EXPERIMENTAL OPEN-LOOP FREQUENCY RESPONSE 
LOAD VALVE SERVO- ACTUATOR 


UNCLASSIFIED 



FREQUENCY, (J , RAD/SEC 

CORRELATION OF CLOSED-LOOP FREQUENCY RESPONSES 
GAS -FARM VALVE AND LOAD VALVE SERVO-ACTUATORS 
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CORRELATION OF OPEN-LOOP FREQUENCY RESPONSES 
TURBINE POWER CONTROL VALVE SERVO-ACTUATOR 


UNCLASSIFIED 


c 

z 

o 

r 

> 

0) 

(/> 

Z3 

m 

□ 


H 


(W 


£ 

►d 


td 

i 

i- 1 


360 r 
340 
320 
300 
280 
260 
240 
220 

<n 200 I 


a iso 

a 


£ 160 
z 

in 140 

in 

< 

^ 120 
100 
80 
60 
40 
20 
0 



CORRELATION OF CLOSED-LOOP FREQUENCY RESPONSES 
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BLEED-LINE VALVE SERVO -ACTUATOR 
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